MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  Of  STANDARDS-  1963-A 


CCAD-H-S3-01 


CENTER  FOR  COMPUTATIONAL  AND  APPLIED  DYNAMICS 
BOSTON  UNITERS  ITT 
BOSTON  MASS.  0221S 


FREE  TAKE  AERODYNAMIC 

ANALYSIS  OF  HELICOPTER  MOTORS 

Laigi  Moriao,  Zaven  Kiptlillu.  Jr.  and  Slobodan  R.  Sipoio 

Mar  «*3 

Final  Report  on 
U.S.  A r*y  Research  Office 
Contract  No.  DAA629-S0-C-0016 

A 

APPROVED  FOR  PUBLIC  RELEASE 
DISTRIBUTION  UNLIMITED 


DTIC 

SELECTE 
JUN  2  1  1983 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  TW»«  PAGE  (BMP  D*lm 

I  REPORT  DOCUMENTATION  PAGE 


12.  «OVT  ACCESSION  NO. 


TR-83-01 

«.  TITLE  fan*  JuKIll.J 

FREE  WAKE  AERODYNAMIC  ANALYSIS  OF 
HELICOPTER  ROTORS 

7-  AUTHOR*-.; 

Luigi  Morino 

Zaven  Kaprielian,  Jr. 

Slobodan  R.  Sipcic 

*■  PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

Boston  University 
Boston,  Mass.  02215 


READ  INSTRUCTIONS 

_ BEFORE  COMPLETING  FORM 

S.  RECIPIENT'S  CATALOG  NUMRER 


S.  TYPE  OP  REPORT  t  PERIOD  COVERED 

FINAL  REPORT  (10/17/79 
_ — _ V?jj/Rl) 

t.  PERPORMINO  ORO.  REPORT  NUMBER 

CC AD-TR- 83-01 

S.  CONTRACT  OR  GRANT  NUMEEfcf*> 

DAAG  29-80-C-0016 


PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  S  WORK  UNIT  NUMBERS 


II.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

v.  t  Army  Rtseareft  Office 
Post  Office  Box  12211 

Research  Triangle  Park,  NC  mtt9 _ 

U.  MONITORING  AGENCY  name  SAOORESSfif  i 


I  IS.  OISTRIBUTION  STATEMENT  (ot  thlm  Mmpmrt) 


12.  REPORT  DATE 

May  26,  1983 

IS.  NUMBER  OP  PAGES 


I  Contntltnt  OlHcm)  I  IS.  SECURITY  CLASS,  (of  thlm  n,ow) 


Unclassified 


IS*.  OECL  ASSI  PIC  ATI ON/ OOWNORAOING 
SCHEDULE 


'Approved  for  public  release :  distribution 
unlimited. 

17.  OISTRIBUTION  STATEMENT  (•!  thm  aA.tr act  nlmj  In  Black  20,  if  dllfmrmt  Ban 


IS.  SUPPLEMENTARY  notes 

INS  VIEW.  OPINIONS.  AND/OR  Ft*IC!NC-3  C0IP7”T0  IN  THIS  REPORT 
ARK  THOSE  OF  THE  AUTHbRIS)  AND  SHCUt  D  C  35  CONSTRUED  AS 
AN  OFFICIAL  DEPARTMENT  OF  THE  ARt/Y  SCD  7‘OM,  POLICY,  OR  DE¬ 
CISION,  UNLESS  SO  DESIGNATED  EY  C7H21  D-'CUMENTATION, 

IS.  KEY  WORDS  (Cmnllmtm  an  nnrM  aMa  if  naaaaaaqr  an*  Hmtllf  *r  Mne*  nuwSarJ 

helicopter,  rotor,  aerodynamics,  wake  dynamics, 
free-wake,  trailing-edge  conditions 


formulation  for  the  free  wake  analysis  of  helicopter 
rotors  in  incompressible  potential  flows  is  presented  here. 
The  formulation  encompasses  both  the  theory  and  its  numeri¬ 
cal  implementation.  For  the  case  of  a  single-bladed  rotor 
in  hover,  the  formulation  is  validated  by  numerical  results 
which  are  in  good  agreement  with  the  generalized  wake  of 
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Landgr ebe.  Extension,  of  the  theory  to  compressible  flows 
is  also  outlined. 

These  results  indicate  that  the  formulation  does  not 
require  any  empirical  assumption  Csuch  as  the  rate  of  con¬ 
traction  of  the  radius  of  the  wake)  in  order  to  avoid  nu¬ 
merical  instabilities.  To  our  knowledge,  the  results  pre¬ 
sented  here  are  the  first  ones  ever  obtained  not  requiring 
any  ad-hoc  assusiptions  in  order  to  avoid  such  problems.  . 
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SECTION  1 


INTRODUCTION 


LuU  Maiiirtirt  Art  j&iifiligti 

A  new  methodology  for  generation  of  the  wake  geometry  for 
the  computational  aerodynamic  aaalyaia  of  a  helicopter  rotor  ia 
preeented  ia  this  report. 

The  availability  of  such  a  methodology  would  enhance 
considerably  the  present  computatioaal  capability  in  this  area. 
This  capability  is  needed  for  instance  for: 

(1)  performance  and  structural  analysis 

(2)  evaluation  of  generalized  forces  for  flutter  analysis 

(3)  evaluation  of  the  outer  potential  velocity  field  for 
the  boundary- layer  and  separated-flow  analysis. 

In  the  classical  rotor-wake  formulation,  the  wake  is 
described  as  a  spiral  (heliooidal)  surface  which  is  obtained  from 
the  assumption  of  uniform  vertical  flow.  This  method  is  not 
sufficiently  accurate  for  the  aerodynamic  analysis  of  helicopter 
rotors.  This  yields  the  need  for  the  development  of  a 
methodology  for  fully-automatic  or  semi-automatic  wake 
generation. 

The  fully  automated  wake  generation  (commonly  referred  to  as 
'free  wake*  analysis)  is  obtained  step-by-step  by  calculating 
from  the  location  of  a  vortex  point  at  a  given  time  step  the 
vocation  at  the  next  timestep:  the  drawback  with  this  approach 
is  that  the  free-wake  analysis  is  quite  expensive  ia  terms  of 
computer  time. 

On  the  other  hand,  a  semi-automatic  wake  generation 
(commonly  referred  to  as  'generalized  wake')  may  be  obtained  by 
expressing  the  analytical  description  of  the  wake  geometry  in 
terns  of  few  parameters  which  are  evaluated  by  fitting 
experimental  results.  The  generalized-wake  analysis  is  accurate 
and  not  more  expensive  than  the  classical-wake  analysis,  but 
currently  requires  the  use  of  expensive  wind  tunnel  experiments 
for  the  generation  of  the  general  ize-wake  model.  On  the  other 
hand,  the  free-wake  analysis  is  more  expensive  than  the 
generalized-wake  analysis  but  less  expensive  than  the  wind  tunnel 
experiments. 

The  objective  of  work  presented  hero  is  the  development  of 
an  efficient  and  general  method  for  free-wake  potential 
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aerodynamic  analysis  which  than  oaa  ba  aaad  (instead  of  the  more 
expensive  experimental  approach)  to  generate  the  general ized— wake 
model  for  nae  in  a  prescribed-wake  analysis.  The  role  of  the 
experimental  work  would  then  be  limited  to  the  validation  of  the 
free-wake  analysis. 


1.2.  Relation  of  Work  with  State  of  the  Art 

An  excellent  review  on  aerodynamic  technology  for  advanced 
rotorcraft  was  presented  by  Landgrebe.  Moffitt.  and  Clark  in 
Refs.  1  and  2.  Additional  reviews  are  presented  in  Refs.  3-7. 
(Compressibility  effects  in  particular  are  reviewed  in  Ref.  7.) 
Therefore  only  works  which  are  particularly  relevant  to  the 
objective  and  the  motivation  of  the  proposed  work  are  included  in 
this  brief  review,  which  is  not  to  be  considered,  by  any  means, 
complete. 

We  feel  that  the  unsteady  analysis  is  the  main  strength  of 
the  method  presented  here  over  other  existing  methods  (accurate 
pressure  evaluation  for  compressible  unsteady  flow  is  badly 
needed  for  flutter  analysis,  see  Ref.  3:  a  correct  and 

efficient  formulation  for  such  a  problem  is  not  available  at 
present  time  even  for  an  isolated  rotor).  Therefore  particular 
emphasis  is  given  in  this  review  to  the  unsteady  flow  analysis. 
Another  important  feature  of  the  methodology  proposed  here  is  the 
capability  of  analyzing  rotor-fuselage  interaction:  this  is  the 
main  advantage  of  so-called  panel  methods  over  lifting-surface 
methods  (see  below).  Therefore  this  aspect  is  also  emphasized  in 
this  review.  A  third  advantage  of  the  methodology  proposed  here 
is  that  it  can  be  extended  to  give  an  exact  integral  equation  for 
nonlinear  compressible  unsteady  flows  with  moving  shock  waves. 
The  extension  of  the  present  formulation  to  compressible  flows  is 
briefly  outlined  in  Appendix  A. 

It  is  interesting  to  examine  the  Concluding  Remarks  made  in 
Reference  2  in  1977: 

'As  with  any  review  of  a  rapidly  changing  subject, 
the  conclusion  that  must  be  drawn  from  this  review  is 
that,  although  much  has  been  achieved  in  terms  of  our 
understanding  the  aerodynamics  of  the  rotor,  much 
remains  to  be  accomplished.  In  recent  years,  there 
have  been  large  advances  in  the  capabilities  of  the 
helicopter  with  today’s  machines  flying  faster, 
farther,  higher,  heavier,  smoother  and  quieter  than 
their  predecessors.  However,  the  gains  possible  with 
today's  advanced  technology  are  becoming  smaller  as  the 
state  of  the  art  approaches  goals  set  using 
yesterday's  ground  rules,  and  the  available  techniques 
must  be  worked  harder  to  achieve  even  modest 
improvements.  Unfortunately,  these  techniques. 
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reflecting  with  a  growing  complexity  an  increased  depth 
of  understanding  of  the  rotor  behavior,  are  beginning 
to  strain  the  capabilities  of  the  competing  facilities 
available.  If  further  advanoea  are  to  be  Bade,  soae 
simplification  (without  loss  of  precision)  is  required. 
This  is  especially  true  if  the  methods  are  to  be  made 
available  to  the  engineer  and  the  designer  at  a 
working,  as  opposed  to  a  research,  level. 

More  specifically,  based  on  this  review,  several 
areas  which  require  more  attention  as  the  industry 
enters  another  decade  of  the  development  of  our 
technology  include:  the  definition  of  the  inflow 
distribution  in  forward  flight  naing  a  generalized  wake 
approach  similar  to  that  developed  in  hover,  the 
combination  of  the  effects  of  unsteady  and  akewed  flow, 
a  resolution  of  the  limits  of  applicability  of  the 
basic  lifting-line  sssuaption,  the  development  of  an 
economic  lifting  surface  analysis,  and  a  representative 
aerodynamic  model  of  the  rotor/airfraae  interaction. 
All  of  these  (and  other  important  elements  not  covered 
in  this  review)  should  be  combined  in  a  tractable,  well 
correlated  analysis,  which  will  be  used  in  combination 
with  an  improved  appreciation  of  the  role  of  the  blade 
aeroelasticity  to  develop  designs  for  tomorrow's 
advanced  rotorcraft.' 


The  work  presented  here  addresses  most  of  these  issues.  As 
mentioned  above,  the  objective  of  this  work  is  to  develop  an 
efficient  and  general  methodology  for  free-wake  potential 
aerodynamic  analysis  of  helicopters  in  hover  or  forward  flight. 
This  methodology  may  be  used  to  generate  a  generalized  wake  which 
is  now  obtained  from  very  expensive  wind  tunnel  results.  (The 
availability  of  such  a  method  and  corresponding  computer  program 
would  enhance  considerably  the  present  computational  capability 
for  an  accurate  evaluation  of  pressure  and  flow  fields  for 
performance  and  structural  analysis.  Such  evaluation  is  also 
required  for  instance  for  the  problem  of  drag-reduction  as  a 
prerequisite  for  the  boundary-layer  analysis.)  Also  mentioned 
above  is  the  fact  that  the  possibility  of  analyzing  compressible 
u n$ta*dy  flows  and  rotor-fuselage  interference  makes  the  integral 
method  presented  here  quite  unique  among  the  'integral  equation' 
methods  whieh  are  well  known  to  be  much  faster  and  much  simpler 
to  use  than  f inite-dif f erenoe  or  finite-element  methods. 

Three  items  which  are  relevant  to  this  report  and  whieh  need 
a  discussion  deeper  than  the  ones  presented  in  Refs.  1  to  7  are 
advanced  computational  methods  (lifting-surface  and  panel 
methods),  wake  roll-up  and  coapressibil ity.  These  items  are 
briefly  examined  in  the  following. 


Consider  first  advanced  computational  methods.  Lifting- 
surface  theories  are  presented  in  Refs.  6  to  7.  Both  methods  are 
of  interest  here.  The  first  (Snmma,  Ref.  6)  introduces  an 
incompressible  time— domain  analysis  (the  wake-dynamics  analysis 
of  this  paper  is  covered  later  in  this  section)  while  the  second 
one  (Sao  and  Sehatzle,  Ref.  7)  introduces  in  a  simplified  form 
(local  Prandtl-Glauert,  chordwise  transformation)  the  effect  of 
compressibility  for  rotors  in  unsteady  flow.  This  paper  is  very 
relevant  here  because  it  clearly  demonstrates  the  importance  of 
an  accurate  evaluation  of  the  wake  geometry  for  the  calculation 
of  the  section  lift  distribution  acting  on  the  blade.  A  third 
lifting-  surface  method  was  developed  by  Suciu.  Preuss  and  Morino 
(Ref.  8)  for  windmill  rotors  and  yields  results  which  are  in 
excellent  agreement  with  those  of  Rao  and  Schatzle  (Ref.  7).  An 
important  development  which  requires  special  attention  is  the 
work  by  Dat  and  Costes  (Refs.  9-12).  who  start  from  the 
acceleration  potential  due  to  a  doublet  and  develop  a  lifting- 
surface  equation  for  compressible  rotor  aerodynamics. 

Next  consider  panel  methods,  a  new  methodology  recently 
introduced  in  aircraft  aerodynamics.  This  methodology  (also 
called  boundary-element  method)  consists  of  the  finite-element 
solution  (over  the  actual  surface  of  the  body)  of  integral 
equations  for  potential  aerodynam ics.  Typically,  the  surface  of 
the  aircraft  is  covered  with  source-panels  (doublet-,  vortex-, 
and  pressure-panels  are  also  used  on  the  surface  of  the  body  and 
of  its  wake).  The  intensity  of  the  source  distribution  is 
obtained  by  assuming  that  the  flow  does  not  penetrate  the  surface 
of  the  body.  (Note  the  difference  with  respect  to  the  lifting- 
surface  formulations,  in  which  the  integral  equation  is  over  the 
mean-surface  of  the  rotor  blades.  By  panel  methods  we  indicate 
only  those  methods  in  which  the  actual  surface  is  used:  their 
main  advantage  over  the  lifting-surface  methods  is  the  capability 
of  analyzing  rot  or- fusel  age  interaction.)  An  early  work  on  the 
flow  field  around  three-dimensional  bodies  by  Hess  and  Smith 
(Ref.  13)  uses  constant  strength  source-elements  to  solve  the 
problem  of  steady  subsonic  flow  around  nonlifting  bodies.  This 
method  has  been  extended  to  lifting  bodies  (Refs.  14-18)  by 
including  doublet,  vortex,  and  lifting-surface  panels.  Work  in 
panel  method  for  unsteady  flow  around  complex  configurations 
include  extensions  of  the  doublet-lattice  method  (Refs.  19  and 
20)  and,  more  recently,  the  work  of  Morino  and  his  collaborators 
(Refs.  21-28). 

This  methodology  has  been  recently  extended  to  helicopter 
aerodynamics.  For  instance,  the  work  by  Dvorak,  Maskew  and 
Woodward  (Ref.  29)  present  a  method  for  calculating  the  complete 
pressure  distribution  on  a  helicopter  fuselage  with  separate  flow 
(the  method  uses  WBAERO  for  the  potential  flow  solution,  the 
boundary  layer  is  calculated  up  to  the  separation  line,  separated 
flow  is  modeled  by  strsamwise  panels  of  uniform  vorticity). 
Additional  results  are  presented  in  Ref.  30.  Straightforward 
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application*  of  paael  aithodi  to  helicopter  aerodynamics  art  also 
presented  by  Soohoo,  Torino.  Noll,  and  Baa  (Safa.  31  aad  32:  a 
formal  a  tioa  for  ooaprassibla  flows  Is  presaatod  ia  Raf.  31). 

Tha  above  remarks  iadieata  that  paa*l-a*rodyaaaies  aathods 
are  baeomiag  available  for  tha  analysts  of  tha  eoaplata 
coaf igaratioa.  Tha  availability  of  saeh  aathods  (aad 
eorraspoadiag  eoapatar  programs)  aahaaeas  considerably  tha 
prasaat  conpata tional  capability  for  an  aeearata  avalaatioa  of 
prassara  aad  flow  fields.  This  avalaatioa  is  baeomiag  more  aad 
more  important  baeaasa  of  raeaat  troads  ia  tha  field  of 
haliooptor  aerodynamics.  Bind  taaaal  experiments  are  vary  costly 
whereas  oompatars  ara  baeomiag  lass  aad  lass  expensive. 
Therefore,  tha  asa  of  oompatars  is  baeomiag  more  attraetiva  for 
tha  aarodyaamie  analysis  of  haliooptar  coaf igaratioa.  For 
iastaaea,  items  saeh  as  higher  performance  (lower  drag,  higher 
spaed,  higher  lift,  higher  reliability)  require  more  theoretical 
analysis:  ia  partiealar,  as  mentioned  above,  redaotioa  of  drag 
requires  a  vary  aeearata  avalaatioa  of  the  potential  flow  field 
as  a  prerequisite  for  tha  bouadary-1 ayer  aad  flow  separation 
analysis  (Raf.  29).  Therefore  panel-aerodynamics  aathods  deserve 
farther  attaatioa  aad  ara  expected  to  beeoma  a  standard  dasiga 
tool  for  airplane  aad  haliooptar  aerodynamics  within  this  deeada. 

Neat  consider  the  issue  of  wake  dyaamies.  An  azeallaat 
review  of  tha  problem  of  tha  wake  roll-up  is  given  ia  Refs.  1  aad 
2  (where  additional  works  aot  iaeladed  her*  are  extensively 
reviewed).  The  essaao*  of  the  state  of  the  art  in  this  area  is 
briefly  summarised  here.  The  various  aerodyaamie  analysis  of  the 
rotor  fall  into  oae  of  the  three  following  types  : 

A.  Classical  wake.  i.e.,  a  wake  geometry  described  by 
a  heliooidal  spiral  with  pass  obtained  from 
uniform  flow  assumption. 

B.  Geaeralised  waks.  i.e.,  a  wake  geometry  obtained 
by  interpolating  experimental  data  ia  terms  of  few 
parameters. 

C.  Free  wake .  i.e.,  s  wake  geometry  obtained 

computationally  as  aa  integral  part  of  the 
solution. 


Aaalytieal  models  for  predicting  the  geometry  of  the  rotor 
wake  were  developed  from  experimental  data  by  Laadgreb*  (Ref.  33 
aad  34),  Crews,  Boheaemser  aad  Ormistoa  (Ref.  35)  aad  Kocurek 
(Raf.  3d).  Laadgrebe's  nodal  was  used  by  Rao  aad  Schatxl*  (Ref. 
7)  ia  their  lifting-surface  theory,  aad  shows  that  a  considerable 
improvement  ia  the  comparison  with  experimental  results  of  Rof. 
37  eaa  be  aad*  simply  by  using  s  geaeralised  wake  geometry 
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instead  of  the  classical  waka  geometry.  Automatic  generation  of 
the  wake  is  considered  for  instance  by  Scully  (Ref.  38).  Summa 
(Ref.  8)  .and  Pouradier  and  Horowitz  (Ref.  39).  All  these  works 
indicate  that  the  algorithms  used  sre  unstable  unless  special 
constraints  (such  ss  specified  contrsetion  ratio)  are  introduced. 
Another  important  issue  is  the  one  of  simplified  algorithms  which 
can  be  used  for  instance  for  modeling  the  far  wake:  several 
models  are  available  for  the  hover  ease  (see  e.g..  Ref.  40-42). 
However  none  of  these  models  is  applicsble  to  the  case  of 
arbitrary  motion  considered  in  this  report. 

As  mentioned  sbove,  it  is  a  generally  accepted  opinion  in 
the  helicopter-aerodynamics  community  that  the  classical-wake 
analysis  is  insufficient  (see  for  instance  Ref.  7)  to  obtain 
accurate  results  and  that  the  free-wake  analysis  is  too  expensive 
and  that,  consequently,  sll  the  effort  should  be  devoted  to  the 
generalized-wake  analysis.  TThile  we  basically  agree  with  this 
assessment,  we  would  like  to  point  out  again  that  in  generating 
the  generalized-wake  geometry  it  is  not  necessary  to  use  the 
costly  wind-tunnel  results:  one  can  use  the  free-wake  analysis 
which  is  more  expensive  than  the  generalized-wake  analysis  but 
also  much  less  expensive  than  the  experimental  approach. 

Next  consider  the  issue  of  compressibility.  As  mentioned 
sbove,  the  analysis  for  compressible  flows  is  examined  in 
Appendix  A.  Hence  for  completeness,  the  issue  of  compressibility 
is  reviewed  here.  The  importance  of  compressibility  was  clearly 
demonstrated  by  Friedman  and  Tuan  (Ref.  3)  for  the  problem  of 
aeroelastic  stability  (i.e.,  flutter  and  divergence)  of  rotor 
blades.  The  work  is  based  on  simple  serodynamic  strip-theories 
(Refs.  43-49).  However,  the  same  effect  is  expected  from  more 
sophisticated  unsteady  three-dimensional  compressible  theories. 
As  mentioned  above,  compressibility  effects  are  included  in  the 
lifting-line  theory  by  Johansson  (Ref.  5),  in  the  lifting-surface 
method  by  Rao  and  Sohatzle  (Ref.  7  which  is  bssed  on  the  work  of 
Rao  and  Jones.  Ref.  50)  and  in  the  works  of  Dat  and  Costes  (Refs. 
9-12)  and  in  the  work  by  Morino  and  Soohoo  (Ref.  31).  A  possible 
alternative  approach  is  the  numerical  solution  of  the 
differential  equation  using  for  instsnce  the  finite-difference 
technique  (Caradonna  et  al..  Refs.  51-54):  this  method  however 
requires  considerable  amounts  of  computer  time.  Therefore  the 
development  of  a  methodology  based  on  the  integral  equation 
presented  in  Appendix  A  appears  to  be  an  attractive  alternative 
to  other  existing  methods  for  the  analysis  of  unsteady 
compressible  flows  around  complex  configurations  (the  vslidstion 
of  this  integral  formulation  is  limited  to  the  particular  case  of 
incompressible  flows:  the  results  are  presented  in  Section  6  of 
this  report). 
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1 .3  .  of  Uork 


The  work  presented  here  include*  two  espect*:  (1)  the 
development  of  *  formulation  for  the  time-dependent  free— wake 
aerodynamic  analysis  of  helieopters  in  hover  end  forward  flight 
and  (2)  validation  of  the  formulation,  for  the  partienlar  eaee  of 
incompressible  flows,  by  developing  a  computer  program  and 
comparing  the  results  against  selected  existing  computational  and 
experimental  data. 

Regarding  the  first  aspect  of  the  effort,  the  formulation 
is  very  general:  the  main  restriction  is  the  assumption  of 
potential  aerodynamics.  This  implies  in  particular  that  viscous 
(attached  and  separated)  flows  are  not  included  here.  For  the 
sake  of  clarity,  compressible  flows  are  dealt  with  in  Appendix  A: 
the  main  body  of  this  report  deals  with  incompressible  flows.  The 
formulation  is  based  upon  sn  integrsl  method  developed  by  Morino 
(Refs.  21  to  23,  see  slso  Section  3  and  S)  for  the  exact 
compressible  three-dimensional  unsteady  velocity-potential 
equation  for  lifting  bodies  hsving  arbitrary  shapes  and  motions. 
(For  incompressible  flows  with  prescribed  wake,  the  formulation 
has  been  applied  to  rotor-fuselage  helicopter  configurations. 
Refs.  31-32,  and  to  time-domain  analysis  of  unsteady  flows  around 
windmill  rotors.  Refs.  8.)  New  theoretical  results  are  included 
in  this  report.  The  formulation  for  wake  dynamics,  presented  in 
Section  2,  is  considerably  different  from  that  of  our  preceding 
publications  (which  did  not  make  explicit  use  of  the  principles 
of  conservation  of  mass  and  momentum  across  the  wake):  we 
believe  it  is  clearer,  more  convincing  and  more  rigorous  than  any 
wake-dynamics  formulation  of  which  we  are  aware.  Also  new  are  the 
results  presented  in  Section  4  which  deals  with  the  issues  of 
wake  generation,  uniqueness  of  solution,  Cotta  condition, 
Joukowski  hypothesis  and  trailing  edge  condition. 

Regarding  the  second  aspect  of  the  effort,  the  validation 
of  the  formulation  includes  time-domsin  free-wake  analysis  of  the 
unsteady  velocity-potential  equation  for  flexible  rotors  in  hover 
or  forward  flight.  The  validation  involves  the  implementation  of 
the  numerical  formulation  into  a  computer  program  and  comparison 
with  existing  experimental  and  numerical  results.  This 
validation  is  limited  to  an  isolated  rotor  in  hover.  However, 
the  formulation  and  the  numerical  algorithm  used  in  the  computer 
program  are  time  aocurate  (i.e.,  they  yield  a  steady  state 
solution  via  an  accurate  time-domain  analysis)  and  therefor*  are 
in  theory  applicable  to  time  dependent  flows  (of  course, 
validation  for  this  application  would  be  required).  The  computer 
algorithm  is  general  in  that  only  the  geometry  and  the  motion  of 
the  surface  of  the  rotor  is  needed  as  sn  input.  The  issue  of  the 
numerical  stability  for  the  automatic  generation  of  the  wake  is 
also  disoussed. 
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SECTION  2 


VAXE  DYNAMICS  IN  INCOMPRESSIBLE  POTENTIAL  FLOWS 


A  general  formulation  for  the  problem  of  the  wake  dynamics 
for  potential  flows  is  presented  in  this  section.  Fundamental 
issoes  related  to  the  wake  dynamics  have  been  clarified  during 
the  effort  reported  here.  These  issues  (such  as  the  boundary 
conditions  on  the  wake  and  at  the  trailing  edge)  are  quite  subtle 
and  have  never  been  discussed  in  our  preceding  publications. 
Therefore  the  theoretical  foundations  of  the  wake  dynamics  are 
presented  here  (see  Sections  2.2  and  2.3).  we  believe  for  the 
first  time.  (Trail ing-edge  issues  are  dealt  within  Section  4.) 
For  the  sake  of  simplicity,  the  formulation  is  presented  for 
incompressible  flows  (compressible  flows  are  dealt  with  in 
Appendix  A). 


lAx.  Incompressible  lay.jtqjj  Flows 

In  this  report  we  will  assume  that  the  frame  of  reference  is 
connected  with  the  undisturbed  air.  We  assume  the  fluid  to  be 
inviscid  and  incompressible.  Hence  the  motion  is  governed  by 
the  Euler  equations  (conservation  of  momentum) 

Dv  1 

~  -  -  -  grad  p  2.1 

Dt  p 

and  the  continuity  equation  for  incompressible  fluid 
(conservation  of  mass) 

div  W  ■  0  2.2 

where  v  is  the  velocity  vector  with  respect  to  a  prescribed  frame 
of  reference,  p  is  the  pressure,  p  is  the  density  (constant  for 
incompressible  fluid),  t  the  time,  whereas 


D_ 

Dt 


3_ 

dt 


+  W*grad 


2.3 


is  the  material  or  substantial  derivative.  Equations  2.1  and  2.2 
form  a  system  of  four  partial  differential  equations  for  four 
unknowns  vx,  vy,  vs,  and  p. 

In  order  to  complete  the  formulation  of  the  problem,  the 
boundary  conditions  at  infinity,  on  the  body  and  on  the  wake  must 
be  obtained. 


r  '  ^ 


Since  the  frsae  of  reference  has  boon  sssuaed  to  bo 
eonnootod  with  tbo  undisturbed  oir.  tbo  boundary  condition  at 
infinity  any  bo  written  aa 

p  ■  pm  and  w  ■  0  (for  P  at  •)  2.4 


On  tbo  body  (rotor  in  oar  esse)  it  ia  aaaaaod  that  tbo 
aurfaoo  of  tbo  body  ia  iaperaeable.  Tbia  iaplioa  that  tbo  noraal 
coapononta  of  tbo  velocity  v  of  tbo  fluid,  and  of  tbo  velocity 
of  tbo  rotor  bladoa  coincide  at  point  P  on  tbo  aurfaoo  ^  of  tbo 
rotor: 


( v  -  ▼fc)*n  •  0  (for  P  on  o^)  2.5 

where  a  ia  tbo  noraal  to  at  P. 

Tbo.  boundary  condition  on  tbo  wake  aro  diaouaaed  in  Suction 
2.4  aftox  introducing  tbo  concopt  of  potential  wako. 


2.2.  Potential  Flows:  Potential  Mil 

Tbo  baaia  of  tbo  diacuaaion  on  tbo  wako  dynaaica  ia  tbo  well 
known  Kelvin's  tbooroa  which  atatea  that  the  circulation 


oror  a  aatorial  contour  C  (i.e.,  a  contour  which  is  aado  up  of 
aatorial  particles)  reaains  constant  ia  tiae.  This  tbooroa  it  an 
iaaodiato  consequence  of  the  definition  of  T.  of  Euler  equations 
(Eq.  2.1)  and  of  the  fact  that  tbo  density  is  constant  (or.  in 
general,  that  the  fluid  is  barotropio). 

Next  asauao  that  tbo  flow  field  is  irrotational  at  tiao  0. 
Then  according  to  Stokes  tbooroa 

1  v.ds  *JJ  curl  v«n  do  2.7 

(whore  C  is  tbo  contour  of  o) ,f~  is  initially  equal  to  aero  for 
any  path  connected  with  a  surface  o  fully  inside  tbo  fluid 
woluae.  Bonce,  for  all  those  paths,  P  reaains  identically  equal 
to  soro.  This  iaplioa  that 


curl  v  ■  0  2.8 

for  alaost  all  tbo  fluid  points  at  all  tiaoa:  tbo  only  points  to 
bo  ozoludod  aro  those  aatorial  points  which  coao  in  contact  with 
tbo  solid  boundaries  (since  for  those  pointa,Kelvin's  tbooroa 
docs  not  apply).  Xa  order  to  siaplify  the  discussion  of  this 
issue,  lot  us  focus  on  the  ease  of  an  isolated  blade  with  a  sharp 


I 

i 
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! 
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trailing  edge  and  consider  only  those  flows  such  that  the  fluid 
leawes  the  surface  of  the  blade  at  the  trailing  edge  (thia  issue 
is  further  discussed  later  in  this  Section  4).  We  call  these 
flows  attached  flows.  Hence  the  points  which  cone  in  contact 
with  the  rotor  are  only  those  emanating  f row  the  trailing  edge 
and  therefore  fora  a  surface:  such  a  surface  is  called  wake. 
Kelvin's  the  ore  a  does  not  apply  to  the  points  of  the  wake.  [Note 
that  if  the  trailing  edge  is  fixed  with  respect  to  the  fraae  of 
reference,  the  wake  is  a  'streak  surface*  (Ref.  52).  Since  the 
trailing  edge  in  general  aowes,  we  say  that  the  wake  is  a 
'generalised  streak  surface.*] 

We  nay  conclude  that  for  an  inwiscid  and  incoapressible 
fluid,  if  a  flow  field  is  initially  irrotatioual.  it  reaains 
irrotational  at  all  tiaes  except  for  those  points  which  coae  in 
contact  with  the  surface  of  the  body.  If  the  flow  is  attached 
(i.e..  by  definition  if  the  fluid  leaves  the  body  only  at  the 
points  of  a  line  which  is  called  'trailing  edge')  then  the  locus 
of  these  points  foras  a  generalized  streak  surface  whieh  is 
called  a  wake. 

It  aay  be  worth  eaphasizing  that  as  shown  above  the  fact 
that  for  attached  flows  the  wake  has  zero  thickness  does  not 
require  futher  assumptions  but  is  a  direct  consequence  of  the 
hypothesis  of  incompressible  (in  general,  barotropic)  invisoid 
flow  and  of  the  definition  of  attached  flows  (the  assumption  if 
initially  irrotational  flows  is  not  essential  for  the  presence  of 
zero-thickness  wakes). 


2  2.  Potent  if 1-Flow  Formulation 

Next  consider  a  well  known  theorem  from  vector  field  theory 
which  states  that  if  a  vector  field  v  is  irrotational  then  there 
exists  a  function,^,  called  velooity  potential  such  that 

V  -  grad  V  2.9 

Hence  our  results  may  be  restated  as  follows:  for  an  inviscid 
incoapressible  fluid,  a  flow  field  whieh  is  attsched  and 
initially  irrotational  is  potential  at  all  tines  and  at  all 
points  except  at  the  point  of  the  wake. 

If  the  flow  is  potential,  i.e.,  if  ▼  is  given  by  Eq.  2.9, 
Eq.  2.1  aay  be  integrated  to  yield  Bernoulli's  theorem 

if  1,  p  1 

—  ♦rlgredpl  -  -p„  2.10 

Ot  2  p  P 


(The  x-derivative  of  Eq.  2.10  yields  the  first  component  of  Eq. 
2.1:  the  constant  on  the  right  hand  aide  of  Eq.  2.10  is  obtained 
froa  the  boundary  conditions  at  infinity.) 
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Furtheraore.  Eq.  22  at;  be  rewritten  as 

V* f  •  0 

vhiti  ▼*  Is  tk«  Laplaciaa  operator. 

Siailarly,  Eq.  2.4  aay  bo  rewritten  as 
(fmQ  (for  P  at  •) 

and  Eq.  2.5  boooaos 


3(£ 

3b 


V 


2J.2 


123 


LL.  Baaa&m  fiMAii&m  aa  Mi 

la  order  to  be  able  to  solve  the  aatbeaatical  problea.  ve 
seed  a  boundary  condition  oa  the  wake.  This  coaditioa  a  ay  be 
obtaiaed  froa  the  priaeiplos  of  conservation  of  aass  aad  aoaeataa 
across  a  surface  of  discontinuity. 

Iadicatiag  by  subscripts  1  aad  2  the  two  sides  of  the  wake, 
let  a  be  the  outward  aoraal  oa  side  1  aad  let 

if  -  fi  -  f,  2.14 

deaote  the  discoatiauity  of  aay  fuactioa  f  across  the  wake 
surface.  (For  the  classical  wiag-vske,  1  aad  2  correspond  to 
upper  aad  lower  sides  respectively,  a  is  the  upper  aoraal  aad  if 

-  -  f  t.) 

The  eqaatioa  of  conservation  of  ass  a  sad  aoaeataa  serose  s 
possible  surface  of  dieooatiauity  (e.g.,  a  shook  wave  or  a  wake) 
are  given  by  (Kef.  58) 


Alp(vn-v#>]  »  0  2J.5 

d[p(vB-vf)v  +  pa]  ■  p  (vB~v#)  Av  ♦  Ap  5  ■  0  2.16 

where  vB  ■  v-n  is  the  aoraal  coapoaent  of  the  velocity  v  whereas 
is  the  velocity  of  the  surface  (by  defiaitioa,  ia  directioa  of 
tha  aoraal  a),  la  the  ease  of  iaeoapressible  flows,  Ap  ■  0  and 
therefore  Eq.  2.15  yields 


Avb  -  0  2.17 

(since  v,  is  the  sane  for  both  sides  of  the  wake).  Then  the 
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normal  component  of  Eq.  2.16  yields 

Ap  «  0  2.18 

Hence  using  Eq.  2.15,  Eq.  2.16  reduces  to 

p(va-v#)Av  ■  0  2.19 

which  implies  either  Av  «  0  (that  is  the  surface  under 
consideration  is  not  a  surface  of  discontinuity)  or  if  there  is  a 
discontinuity  (wake  surface) 


vn  -  v,  2.20 

that  is  the  flow  does  not  penetrate  the  surface  of  the  wake. 

Note  that  Eq.  2.20  implies  in  particular  that  only  the 
tangential  components  of  the  velocity  are  discontinuous  across 
the  surface  of  the  wake,  that  is  (using  a  local  frame  of 
reference  ?a,  T*,  T*  with  T,  in  the  direction  of  the  normal  n) 


Av  ■  Aivjij  +  v*!*  +  v_n) 


A(i*  ~ix  +  it) 
9xx  9 i| 


Jlf  ix  *  T*  “  »rBdtAf  2J1 


where 


9t  -  d f  - 

gradtf  •  r~  ii  +a— -  1*  2.22 

1  011  OX| 

is  the  'tangentisl  portion'  of  the  vector  grad  f.  (It  is 
important  to  note  that  grad  Af  is  a  meaningless  expression  since 
Af  is  defined  only  over  the  surface  of  the  wake.) 

Next  consider  Bernoulli's  theorem  Eq.  2.10  with  the  wake 
condition.  Ap  »  0,  Eq.  2.18.  This  yields 

3«£i  _  £f£i  +  i(va.  vi  -  ?*-v»)  •  0  2.23 

at  at  2 
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that  is 


<f»)  +  -(Vi  +  Vi)‘(vx  - 


V») 


*  2' l5*  *  ■ 0 


2.24 


or,  introducing  tho  velocity  of  a  wake  point,  Pw, 


%  -T^‘  +  *»> 


225 


and  using  Eq.  2.17, 


(—  ♦  vw-gradt)At^  -  0 


2.26 


The  above  equation  ia  the  desired  boundary  condition  on  the 

wake,  and  nay  be  rewritten  as 


Ot 


•A«p  *  0 


2.27 


where 


Dw  a  a  a  a 

—  ■  r-  +  v’grad,  -  —  ♦  v„.  —  +  v_  - — 
Dt  3t  w  1  3t  ,l  3*i  w*dx. 


2.28 


is  the  aaterial  tine  derivative  for  a  function  defined  only  over 
a  surface:  in  a  frsae  of  reference  aoving  with  (constant) 
velocity  equal  to  vy  at  a  given  tine  and  a  given  point,  the 
substantial  derivative  coincides  with  the  partial  tine 
derivative.  This  iaplies  that 

•  constant  229 

following  a  point  P_  which  has  velocity  t.  given  by  Equation 
2.23. 


Note  that  v  and  are  the  average  values  of  the 

coaponents  1  and  2  respectively  of  vt  and  v t  (which  are 
discontinuous  across  the  wake),  whereas  vva  ■  vn  ■  v,  (see 
Equation  2.16)  since  v&  is  continuous  across  the  wake.  The  above 
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results  can  b«  rtititid  ia  aon  physical  tataa  saying  that  a 
poiat  of  tha  wa ha  has  a  velocity  v  equal  to  tha  average  batvaaa 
tha  velocities,  vx  a ad  5,,  oa  sidas  1  sad  2  of  tha  vaka 
respectively.  Tha  value  of  Af(Pw)  is  tiaa  iadapaadaat  sad  equal 
to  tha  value  that  it  had  whea  tha  poiat  Pv  laft  tha  trailing 
adga. 


It  nay  ba  worth  aotiag  that  wva  doos  aot  appear  ia  Equation 
2.24.  Bowawar  tha  fact  that  Pw  has  aoraal  velocity  aqaal  to  twb 
ataas  froa  tha  fact  that  tha  poiats  of  tha  surface  hare  Telocity 
aqaal  to  t,  »  tm»  t,.  (Eqaatioa  2 .20) .  This  is  ia  agreement 
with  tha  results  of  Section  2.2  that  tha  wake  surface  is  a 
gaaarslisad  streak  surface. 


1.3.  Summary 

la  summary,  if  tha  flow  field  of  aa  inwiseid  incompressible 
fluid  is  initially  irrotatioaal.  it  remains  irrotatioual  at  all 

times  sad  therefore  may  ba  given  aa 

v  -  grad  2J0 

Than  (sea  2.11),  f  satisfies  Laplace  equation 

▼*<f  -  0  2.31 

with  coaditioa  at  infinity,  (Eq.  2.12) 

Cf  •  0  (for  P  at  •)  2.32 

Oa  the  other  head  (see  Bq.  2.13),  the  bouadary  condition  on 
the  rotor  blade  is  given  by 


HM  V“  2J3 

whereas  oa  the  wake  (see  Eq.  2.29) 

dy<Pw)  ■  constant  ia  time  2.34 

Equations  2.30  to  2.34,  with  the  addition  of  Joukowski's 
hypothesis  (see  Section  4),  may  ba  used  to  obtain  tha  solution 
for  <f.  Oaea  <f  is  kaowa,  the  perturbation  velocity  may  ba 
avslustad  using  Eq.  2.30.  Than  tha  pressure  may  ba  avalustsd 
using  Bernoulli's  theorem 


2.35 
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SECTION  3 


INTEGRAL  FORMULATION 


The  integral  equation  used  in  this  work  is  a  particular 
case  of  that  introduced  by  Morino  (Refs.  21-23)  for  the  general 
case  of  potential  compressible  flows  for  bodies  having  arbitrary 
shapes  and  motions.  For  the  sake  of  completeness  and  clarity  the 
derivation  of  snch  an  equation  is  briefly  outlined  here  for  the 
specific  case  of  interest*  incompressible  potential  flow  around  a 
helicopter  rotor  having  arbitrary  shape  and  motion.  The  integral 
equation  is  based  on  the  classical  Green's  function  method.  This 
requires  the  use  of  Green's  formula  and  Green's  theorem  which  are 
presented  first. 


3.1  Green's  Formula  for  Laplace's  Equation 

Consider  two  arbitrary  functions  f  and  g  and  note  that 

div(f  gradg)  ■  gradf -gradg  +  f  V2g  3.1 

div(g  gradf )  -  gradg -gradf  +  g  V*f  3.2 

div(f  gradg  -  g  gradf)  *  f  V*g  -  g  V*f  3.3 

Also,  according  to  Gauss  theorem,  for  any  arbitrary  vector  a 

JIfv,  div  1  dV  "  -fyo  5'5  da  3-4 


and 

or 


where  a  surrounds  V  and  n  is  the  inward  directed  normal. 
Combining  Eqs.  3.3  and  3.4  yields 

”  JJ*J  (f  V  g  -  g  V  f)dV  *  (f  gradg  -  g  gradf)- n  do 


#  “  H  -  £ 


3.5 


Next,  assume  that  both  f  and  g  satisfy  the  Laplace  equation, 
that  is 


V*f  -  V*g  -  0  (in  V). 


3.6 


3-1 


Then  Equation  3.5  yields 


jj. (f  &  -•  H> d*  ■ 0  3-7 

which  is  the  desired  Green's  formula  for  the  Laplace  equation. 

Next  choose  the  function  g  to  be  the  simplest  function 
which  satisfies  the  Laplace  equation  with  the  condition 

g  =  0  at  infinity  3.8 

If  g  is  a  function  of  the  distance  r  from  a  specified  point  P*» 
then  the  most  general  function  that  satisfies  the  Laplace 
equation  is 

g  *  A/r  +  B  3.9 

where 

r  =  IP  -  P*|  3.10 

Using  Equation  3.8  yields  B  =  0.  Note  that  g  3  A/r  represents  the 
velocity  potential  of  a  source.  It  is  convenient  to  choose  A 
such  that  g  represents  a  unit  source,  that  is,  a  source  with  flux 
equal  to  one.  The  flux  through  a  spherical  surface  of  radius  R 
is  given  by 

v.n  do  3  tP  — 1  do  =  -  4nR*  3  1  3.11 

ar|r»R  R 

which  implies  A  3  -l/4n  and 


8 


1 

4nr 


3.12 


Combining  Eq.  3.7  and  3.8  yields  the  desired  Green  formula 
for  the  Laplace  equation 

ff [f  r-(r~)  '  7T  (r~)]  d*  “ 0  3-13 

an  4nr  on  4nr 

Note  that  P,  must  be  outside  the 
is  not  satisfied  at  P,. 


volume  V,  otherwise  Equation  3.6 


3-2-  green's  Theorem  Xfi*  Incompressible  Potential  Aergdnaaiat 

In  order  to  obtain  Green's  theorem  for  incompressible 
aerodynamics,  let  the  function  f  in  Eq.  3.13  be  the  velocity 
potential  .  Hence  the  volume  V  must  exclude  the  volume  of 
the  rotor  blades  as  veil  as  a  thin  layer,  Vv,  which  includes  the 
vake  surface  ov  (because  the  equation  V  <p  *  0  is  not  valid  in 
and  on  ow).  In  addition,  the  volume  V  must  exclude  a  small 
volume  around  P.  because  the  equation  T  |  *  0  is  not  valid  at  P,. 
Finally,  the  volume  V  must  be  bound.  Hence  let  the  surface  a  be 
the  boundary  of  the  volume  V.  This  surface  is  the  sum  of  three 
surfaces:  the  surface  o^w  (which  surrounds  the  rotor  blade 

volume,  Vj,,  and  the  wake  volume,  Vw.*  note  that  the  inward  normal 
for  a  is  outward  to  o^),  the  surface  oe  (which  is  a  spherical 
surface  of  radius  s  and  center  P,)  and  the  surface  om  (which  is  a 
spherical  surface  of  radius  R  and  center  P*). 

Then  Eq.  3.13,  for  £  m  if.  reduces  to 


d<p  -1  d  -1 

[  If  (-_-)  -  if  —(-—)] 

.  dn  4nr  3n  4«r 

Note  that,  as  the  radius  s  of  o,  goes  to  zero. 


3.14 


ft 


1  is 
*♦0  JJ<xm 


(  (li->  -  if  ~(~)]  do 

a-  A—  1  3n  4nr 


llf’  spft.1'  *  jnft.  k  J<’1 


“‘f* 

whereas,  as  the  radius  R  of  o„  goes  to  infinity 


lim 


3  tf 

3n 


)]  do 


~  lim  (  if  + 
R 


(i^j  Ri 

3n 


0 


under  the  condition 


lim  f  •  0 
R-® 


3.1,5 


3.16 


3.17 
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3.18 


lia  ^  E  ■  0 
3a 


[It  is  verified  a  posteriori  (Eqs.  3.35-3 .36)  that  these 
conditions  are  satisfied.] 


Therefore,  as  e  tends  to  zero  and  E  tends  to  infinity,  Eq. 
3.14  becoaMS 


3.19 


It  nay  be  worth  noting  that  if  the  point  P*  is  inside  a.,, 
then  there  is  no  need  for  the  surface  og.  In  this  case.  Eq.  3.14 
yields  siaply 


0.#  (r^-)  “  ^ — >  1  do 

JT-  da  4itx  1  da  4  nr 
®bw 

Equations  3.19  and  3.20  nay  be  coabined  by  writing 

E-r-  -  # 


bw 


where 


E,  *  E(P«)  *  1  (P*  outside  o^) 

■  0  (P*  inside  o^) 


3.20 


3.21 


3.22 


3.3.  IlJu.  Contribution 

Next  let  ov  becoae  infinitesimally  close  to  the  surface  of 
the  wake.  Note  that  in  this  process,  the  closed  surface,  ow, 
the  wake  is  replaced  by  the  two  sides  of  an  open  surface  o^.  Let 
n  be  the  normal  on  the  side  1  of  ov  .  In  the  liait  one  obtains 

£ ; - '  If,;  "S;  *  ‘ 0 

[since  (see  Eq.  2.20)  A(^)  *  Avn  »  0]  whereas 
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3.24 


ft.*  b(;>  ‘fj;';1  d° 

°w  w 

Therefore,  in  the  limit.  Eq.  3.21  reduces  to 

4«E*f.  -  -#  $  7  -  f  k(;»  d»  *11-  **  n';’ dd  3-25 

o^,  ow 

» 

where  o^  is  the  (closed)  surface  of  the  rotor  blade  and  ow  is  the 
(open)  surface  of  the  wake  of  the  rotor  blade.  Furthermore  (see 
Eq.  2.14).  A  if*  ‘fi  ~  <f*  .  whereas  n  is  the  normal  on  the  side  1  of 
the  wake.  Note  that  A <f  on  <x  is  evaluated  from  Eq.  2.34.  Note 
also  that  the  vortex-layer  waxe  of  the  rotor  is  represented  as  a 
doublet  layer.  The  proof  of  the  equivalence  of  doublet  layers 
and  vortex  layers  is  given,  for  instance,  in  Reference  57. 


3 .4 .  Integral  Equation 

Equation  3.25  may  be  used  to  obtain  the  value  of  %f  at  any 
point  in  the  field  if  the  value  of  cf  and  d(f/ dn  on  the  surface 
of  the  rotor  and  the  value  of  A<f  on  its  wake  are  known.  Note 
that  df/dn  is  known  from  the  boundary  condition  on  the  rotor,  Eq. 
2.33,  whereas  Af  may  be  calculated  from  the  boundary  condition 
on  the  surface  of  the  wake,  Eq.  2.34. 

Hence,  in  order  to  be  able  to  use  Eq.  3.25  one  must  have  an 
equation  to  evaluate^  on  the  surface:  such  an  equation  is 
obtained  by  noting  that  if  P,  approaches  a  point  on  the  surface 
of  the  rotor,  then  the  value  of  approaches  the  value  of  « f  on 
the  surface  of  the  rotor.  In  order  to  perform  this  limit,  it  is 
convenient  to  interpret  the  doublet  integral  in  terms  of  solid 
angles.  Note  that  (see  Figure  3.1) 


In  particular,  for  a  closed  surface  and  1,  Eq.  3.26  becomes 
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&  r(r”)  do  ■  iD  -  Q* 

Jjg  an  4nr  AitJJQ  4ji  * 

Note  that,  as  is  wall  known 

Q*  *  1  for  P*  inside  o^ 

■  0  for  P*  outside  o^ 

Coopering  Equation  3.22  and  3.28  yields 

E«  *  1  -  Q*/4n  (P,  inside  and  ontside  a^>  3.29 

Using  Equations  3.29,  Equation  3.25  may  be  rewritten  as 

4n<f*  -  -£j>  ~*i  (i)  do 

■''or.  3n  r 
D 

+  (  f -  <f*)n<  r)dff  +  Af  rn(ir)do  3-30 

The  advantage  of  Equation  3.30  over  Equation  3.25  is  that 
each  tern  is  continuous  (as  P*  crosses  or^)  and  therefore  is  valid 
in  particular  even  if  P,  is  on  o^. 

Having  established  this,  one  can  now  go  back  to  Equation 
3.25  with  E«  given  by  Equation  3.29  which  is  valid  also  for  P*  on 
0^  (and  coincides  with  Equation  3.22  if  P*  is  not  on  o^).  In 
particular  if  P*  is  a  regular  point  of  o^  (i.e.,  a  point  where  ob 
has  a  unique  tangent  plane  )  then  ■  2n  and  E«  *  1/2.  Equation 
3.22  nay  thus  be  generalized  as 

E*  *  1  -  Qm/2v  *1  P,  outside 

■  1/2  P,  on  o^  (regular  point) 

"0  P*  inside  3.31 

In  any  event,  for  P*  on  Ov,  Equation  3.25  (with  E»  given  by 
Equation  3.29,  or  Equation  3.il  if  P*  is  a  regular  point  of  o^) 
is  an  integral  equation  relating  the  unknown  values  of  the 
velocity  potential  on  the  surface  of  the  rotor,  to  the  values  of 
d<f/dn  (prescribed  by  the  boundary  condition  on  the  surface  of  the 
blade)  and  the  values  of  the  potential  discontinuity  on  the  wake 
(known  from  the  preceding  time  history). 

The  time  evolution  described  by  Equation  3.25  and  its 
interpretation  in  physical  terms  is  presented  in  Section  4  along 
with  a  discussion  of  the  trailing  edge  condition. 


3.27 


3.28 
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Li  Conditions  H  jgliaiJg 

la  this  sabsaetioa  Equations  3.1?  sad  3.18  ars  vsrifiad. 
Not#  as  r  goas  to  infinity.  Equation  3.43  yialds 


*■  #.  ‘sf '  jb>  - f 


*bt 


'  th'ftjir  a°' 

®bu 


*bw 


Sine*  tha  flu  through  o^w  is  aqul  to  zaro,  i.a. 


3.32 


(this  condition  oust  ba  satisfiad  by 
using  g  ■  1  in  Equation  3.7),  wharaas 


u 

da 


3.33 


and  is  obtainad  by 


lia  ff  do  ■  fiaita,  3.34 

r 


Equation  3.32  iaplias 

■  0(r"*)  at  • 


3.33 


and 


a 

an 


Ofr**)  at  • 


3.36 
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SECTION  4 


WAKE  GENERATION 


la  this  section  the  issne  of  the  vake  generation  is 
analyzed.  First  the  analytical  results  of  Section  3  are 
interpreted  fros>  a  physical  point  of  view.  Next  Joukowski 
hypothesis  (of  smooth  flow  at  the  trailing  edge)  is  introduced 
and  then  the  wake  generation  and  the  trailing  edge  condition  are 
discusaed. 

We  believe  in  particnlar  that  the  issne  of  the  trailing 
edge  condition  (i.e.,  Kutta  condition,  Joukowski  hypothesis)  is 
not  well  understood  and  quite  confusing  in  the  existing 
literature  (and  it  may  be  a  source  of  error  in  the  numerical 
formulation).  Whereas  this  issue  is  discussed  in  greater  details 
in  Reference  S3,  the  major  conclusions  reached  there  are 
summarized  here  in  the  hope  of  shedding  some  light  on  this  aspect 
of  the  foundations  of  potential  flows. 


ixi  Sudden 

In  order  to  discuss  the  problem  of  wake  generation  it  is 
convenient  to  consider  the  problem  of  a  rotor  subject  to  sudden 
start,  that  is  the  case  of  a  rotor  which  for  t  <  0  is  at  rest 
and  is  surrounded  by  a  fluid  which  is  also  at  rest.  At  time  0* 
the  rotor  is  subject  to  a  motion  (which  for  simplicity  we  will 
assume  to  be  of  rigid-body  type)  with  finite  velocity. 


Since  for  time  t  <  0  the  fluid  was  at  rest,  the  wake  has 
not  yet  had  time  to  form.  Therefore  ■  0  everywhere  in  the 
field  and  the  solution  to  the  problem  at  time  0*  is  obtained  by 
solving  Equation  3.2S  with  “  0,  that  is 


4,E.f.  .-$>  i|*;-  f  -<;»  ■»« 

ffb 


4.1 


with  E*  *  1  -  Q*/4n  (E*  »  1/2  on  regular  points  of  o^,  see 
Fq.  3.31).  This  corresponds  to  the  solution  of  Laplace  equation 
without  wake  contribution  and,  as  well  known  (for  instance  from 
conformal  mapping  solution  in  two  dimensional  flows),  this  yields 
a  separation  point  which  is  different  from  the  trailing  edge  and 
hence  infinite  velocity  at  the  trailing  edge.  In  turn  this 
indicates  the  presence  of  a  vortex  at  the  trailing  which  is 
manifested  by  the  fact  that 


Afte  *  0 


4.2 
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{This  can  be 


seen  from  the  fact  that  according  to  Stokes  theorem 


where  C  is  a 
intensity  of 


-  r 


contour  around  the  trailing  edge  whereas 
the  vortex  at  the  trailing  edge.] 


4.3 


f  is  the 


4.2  Uniqueness  of  Solution:  Kutta  Condition 

Perhaps  less  well  known  is  the  fact  that  the  solution  for 
the  problem  is  unique  (the  proof  of  uniqueness  is  given  for 
instance  in  Batchelor,  Ref. 57).  More  precisely  it  may  be  shown 
that  if  the  solution  to  exterior  Neumann  problem  for  Laplace 
equation  is  single-valued  then  it  is  also  unique.  Multi-valued 
solutions  cannot  occur  in  simply  connected  regions.  This  is  in 
sharp  conflict  with  the  classical  two  dimensional  results  in 
which  an  arbitrary  vortex  can  be  added  inside  an  airfoil  to 
obtain  smooth  flow  at  the  trailing  edge.  The  mathematical  reason 
for  the  difference  is  that  the  flow  region  around  an  airfoil  is 
doubly  connected,  whereas  that  around  a  rotor  (as  well  as  most 
three  dimensional  shapes  of  aerodynamic  interest)  is  simply 
connected.  From  a  physical  point  of  view  we  can  say  that  it  is 
possible  to  add  a  vortex  of  arbitrary  intensity  inside  an  airfoil 
(or  inside  a  doughnut-shaped  object,  which  also  yields  a  doubly 
connected  flow-region):  then,  adding  a  suitable  single-valued 
solution  one  obtains  a  nontrivial  solution  to  the  homogeneous 
Neumann  exterior  problem  for  Laplace  equation.  This  solution  can 
always  be  added  to  the  solution  of  the  airfoil  problem,  which  is 
therefore  nonunique.  (From  a  mathematical  point  of  view  it  may  be 
noted  that  the  field  caused  by  the  vortex  is  multivalued  and  the 
proof  of  uniqueness  fails  for  multivalued  potential  functions). 

In  the  case  of  a  rotor  it  is  possible  to  add  a  vortex  inside 
the  rotor-surface:  however  it  is  not  possible  to  have  a  contour 
which  is  'interlocked*  with  the  vortex  without  penetrating  the 
rotor  surface  (this  is  true,  by  definition,  for  any  simply 
connected  region).  The  potential  flow  is  then  single-valued  and 
is  impossible  to  generate  a  nontrivial  solution.  The  solution  to 
the  rotor  problem  is  thus  unique. 

In  order  to  clarify  the  issue  of  the  trailing  edge 
conditions  we  will  call  Kutta  condition  as  that  trailing  edge 
condition  which  is  used  to  eliminate  the  issue  of  a  nonunique 
solution.  This  condition  is  required  for  two-dimensional  flows, 
but  not  in  three-dimensional  flows  (unless  the  region  is  multiply 
connected)  in  which  the  wake  is  responsible  for  the  elimination 
of  singularities  at  the  trailing  edge.  From  a  mathematical  point 
of  view  the  Kutta  condition  requires  that  the  vorticity 
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distribution  over  the  airfoil  dots  not  go  to  infinity  at  the 
trailing  adga.  Th i s  is  different  from  the  condition  that 
discrete  vortices  cannot  remain  at  the  trailing  edge:  each  a 
condition  is  a  consequence  of  Joukowski  hypothesis  discussed  in 
Section  4.3. 


±xl  Joukowski  glB&itilit 

The  asansiption  introduced  here  is  that  the  phyeical  behavior 
of  the  fluid  is  such  that  the  'flow  is  smooth  at  the  trailing 
edge.'  Such  an  assumption  is  connected  with  the  works  of 
Joukowski  (Ref.  59)  and  Eutts  (Ref.  60).  In  order  to  avoid 
confusion  with  the  Eutts  condition  discussed  in  Section  4.2.  we 
will  refer  to  the  above  assumption  as  Joukowski  hypothesis  (for  a 
historical  account  on  the  works  of  Eutta  sad  Joukowski  see  Refs. 
57  and  61). 

In  particular  we  assume  that  if  concentrated  vortices  form  at 
the  trailing  edge  (this  may  happen  because  of  a  sudden  start,  as 
discussed  in  Section  4.1,  as  well  as  any  other  time-discontinuity 
in  the  velocity  distribution  on  the  surface  of  the  rotor),  then 
these  vortices  are  immediately  shed.  At  any  other  time  there  are 
no  concentrated  vortices  at  the  trailing  edge:  the  implication  of 
this  assumption  is  that  the  value  of  is  continuous  at  the 
trailing  edge,  i.e.,  if  Pw  is  a  point  on  the  wake  and  Pte  is  a 
point  on  the  trailing  edge, 

lim  A<f(Pw)  *  4.4 

Pw-Pt. 

where  A<pt#  is  the  difference  of  the  values  of  the  potential  as 
the  points  on  upper  and  lower  two  sides  of  the  rotor  blade 
approach  a  point  of  the  trailing  edge. 


±A  Trailing  Edge  Plow 

In  order  to  appreciate  the  Joukowski  hypothesis  as  it  relates 
to  the  issue  of  uniqueness  of  the  solution,  let  us  examine  in 
detail  the  flow  at  the  trailing  edge.  In  particular  we  want  to 
examine  the  limiting  case  of  the  wake  boundary  condition,  Eq. 
2.16,  as  a  wake  point  approaches  the  trailing  edge.  This 
analysis  is  an  extension  to  unsteady  flows  of  the  work  by  Mangier 
and  Smith,  Ref.  62,  but  the  interpretation  of  the  results  is 
quite  different  from  theirs.  Let  be  a  unit  vector  tangent  to 
the  trailing  edge  end  ev  be  a  unit  vector  in  the  plane  of  the 
wake  and  normal  to  ex(«s9  Figure  4.1).  Let  vT  and  vv  be  the 
components  of  vw  ia  direction  et  and  e^ 
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4.5 


The  sense  of  e  is 
the  rotor,  so  that 


-  J  <vt  +  v,)  -  vrer  +  iJ 
arbitrary,  bnt  the  sense  of  e^  is  away  from 


v 


*> 


> 


0 


Also  let  &r  and  6^  be  the  components  of 

Av  -  vx  -  v,  «  6Tet  +  6^  ey 


4.6 


4.7 


With  these  definitions,  Eq.  2.26  may  be  written  as 

A<f  +  vt5t  +  v,8^  ■  0  4.8 

It  may  be  worth  noting  that  Eq.  4.6  is  the  limiting  case  of 
Eq.  2.26  which  may  be  integrated  to  yield  Eq.  2.29,  which 
describes  the  fact  that  Aq>  is  convected  by  the  flow.  Therefore 
we  believe  that  Eq.  4.6  is  not  an  additional  condition  to  be 
imposed  but  is  a  relationship  that  is  satisfied  since  the 
solution  satisfies  Eq.  2.29.  However  such  a  relationship  has 
interesting  properties  that  help  clarify  the  issue  of  Joukowski 
hypothesis. 

The  case  of  interest  here  is  that  of  a  rotor  with  a  finite 
trailing-edge  angle.  Considering  Figure  4.2,  there  are  five 
cases  of  interest  here.  In  the  first  two  cases  (a  and  b)  the 
velocity  is  infinite  at  the  trailing  edge:  according  to  Joukowski 
hypothesis,  a  concentrated  vortex,  located  at  the  trailing  edge, 
is  immediately  shed.  Therefore  this  case  occurs  only  at  a 
discrete  number  of  times  and  thus  the  other  three  cases  are  the 
only  ones  which  need  further  discussion.  Consider  first  the  case 
(c):  the  velocity  on  both  sides  of  the  rotor  must  have  zero 
components  at  the  trailing  edge  (in  a  frame  of  reference  having 
the  velocity  of  the  trailing  edge  point  0)  along  the  normal  to 
the  rotor,  because  of  the  body  boundary  conditions,  Eq.  2.5  (this 
is  obtained  automatically  with  the  integral  representation  given 
by  Eq.  2.25)  and  must  also  have  zero  components  (in  that  same 
frame  of  reference)  in  direction  normal  to  the  wake  (which  is  a 
generalized  streak  surface).  If  two  components  are  equal  to  zero 
then  the  whole  vector  of  the  cross-flow  velocity  (velocity  of  the 
flow  yield  in  the  plane  of  the  figure,  i.e.,  normal  to  the 
trailing  edge)  is  equal  to  zero.  Hence 

0  4.9 

whereas  v,^>  0.  Note  that  according  to  Eq.  4.4, 
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4.10 


+  ▼»„>  -  \  v,y  >  0 

whereas 

5y  *  vXy  ■  V*V  *  "v»y  <  0  4.11 

Therefore  cate  (c)  can  occur  if  and  only  if 

A^  +  ▼rsT  >  0 

Similarly,  case  (d)  implies  vt)>  >  0,  v,v  *  0,  and  hence 

▼w  -  jvly>  0 

SV  “  vty>  0 

and  therefore  such  case  is  possible  if  and  only  if 

A<f  +  vTSt  <0  4.14 

Finally  consider  case  (e):  in  this  case  v4  and  v*  are  both 
zero.  Hence  both  vv  and  5^  are  equal  to  zero  and  Eq.  4.6  reduces 
to 

+  vt5t  »  0  4.15 

The  flov  moves  parallel  to  the  trailing  edge,  and  the  vorticity 
is  perpendicular  to  the  trailing  edge.  [Note  that  for  two- 
dimensional  flows  Vj  «  &t«  0  which  implies  A«f«  0:  hence  for 
two-dimensional  flows,  case  (e)  occurs  if  and  only  if  A«j»  0 
(e.g..  steady  state).  Case  (c)  occurs  if  A<f>  0  and  case  (d) 
occurs  if  Af<  0.] 

LlL  Pfcytltal  Inte rare tat  ion 

Let  us  go  back  to  the  case  of  impulsive  start.  A  vortex  is 
generated  at  the  trailing  edge  and  immediately  shed.  Then  the 
wake  is  continuously  shed  from  the  trailing  edge  and  convected 
with  velocity  v^  *  (v*  +  v»)/2.  The  wake  is  typically  tangent  to 
one  of  the  two  sides  of  the  rotor.  This  causes  a  cross-flow 
stagnation  point  (high  pressure)  on  one  side  and  a  low  pressure 
point  on  the  other  side. 

It  nay  be  worth  emphasizing  that  in  deriving  the  above 
formulation  we  have  nsed  only  the  principles  of  conservation  of 
mass  and  momentum,  plus  the  assumptions  of  inviscid, 
incompressible  and  initially  irrotational  flows  and  the  Joukowski 
hypothesis.  The  correct  'amount  and  distribution’  of  vorticity 
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4.12 


4.13 


(as  doublet  layers)  is  predicted  by  these  basic  assumptions.  It 
may  be  said  as  the  viscosity  goes  to  zero  the  solution  to  the 
Navier  Stokes  equations  is  that  of  Euler  equation  with  the 
addition  of  Joukovski  hypothesis.  (As  mentioned  above,  Kutta 
condition  of  finite  vorticity  distribution  at  the  trailing  edge 
is  needed  to  insure  the  uniqueness  of  the  solution  only  for  two- 
dimensional  flows:  the  nonuniquene ss  arises  only  for  multiply- 
connected  regions,  as  the  one  of  a  two-dimensional  airfoil.) 
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SECTION  5 


COMPUTATIONAL  ALGORITHM 

The  computational  algorithm  used  to  obtain  an  approximate 
solution  of  Eq.  3.25  is  presented  in  this  section.  The  surfaces 
of  the  rotor  and  of  its  wake  are  divided  into  small  surface 
elements,  <j..  The  potential  and  the  normalwash  are  assumed  to  be 
constant  with  each  element.  This  yields  a  set  of  equations 
relating  the  values  of  the  potential  if.  at  the  centroids  of  the 
elements,  a.,  to  the  values  of  the  normalwash  at  the  centroids  of 
the  elements,  a-  and  to  the  values  of  at  the  centroids  of  the 
wake  elements.  .  The  wake  geometry  is  evaluated  step-by-step 
as  follows:  the  velocity  of  the  fluid  at  the  corners  of  each 
element  is  evaluated  and  from  this  the  location  of  the  element  at 
the  new  time  step  is  obtained. 


5 .1 .  Theoretical  Algorithm 

The  results  of  the  preceding  sections  may  be  summarized  as 
follows:  assume  that  at  time  t  =  t,  the  wake  geometry  and  the 
potential-discontinuity  distributions  are  known,  then  the  value 
of  the  potential  is  obtained  by  solving  Eq.  3.25 

E-f-  '  If  d" 

-  ||  Aif  — (— )  da  5.1 

JJa,  ~  da  4nr 
w 

with 

4 >  =  dJ£  5.2 

T  3n 


and  E*  given  by  Eq.  3.29.  Once  is  known,  Eq.  5.1  may  be  used 
to  calculate  the  veloci./  at  any  point  in  the  fields  as  v,  = 
grad,  ,  (where  grad*  indicates  differentiation  with  respect  to 
P»).  Noting  that 

grad,(-)  =  -,  5.3 
r  r 


where 


r  =  P  -  P* 


5.4 
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one  obtains 


-<f£  tt~*  "  ^-(-^-,)]do 

JJa  4nr  '  da  4nr 


JL  j^7‘ > do 
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This  equation  may  also  be  used  to  calculate  the  velocity  of  the 
points  of  the  wake  if  the  contribution  of  the  wake  integral  in  an 
infinitesimal  neighborhood  of  the  wake  point  is  excluded  from  the 
calculation  (such  a  contribution  is  responsible  for  the  velocity 
discontinuity  across  the  wake  and  its  exclusion  automatically 
yields  the  semi-sum  between  the  two  values  on  the  two  sides  of 
the  wakes). 

Once  the  velocity  vw  (of  a  wake  point  Pw)  is  evaluated,  the 
new  position  of  the  wake-point  Pw  may  be  evaluated  as 


Pw(t,+  dt)  =  P„(t,)  +  vw(Pw,t*)  dt 
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Note  that  a  new  infinitesimal  strip  of  wake  surface  has  been 
added  (between  the  trailing  edge  and  the  locus  of  the  points 
which  were  on  the  trailing  edge  at  time  t*).  Note  also  that, 
according  to  Eq.  2.34,  A>f(Pw)  does  not  change  in  time  (i.e..  Ay 
is  'convected'  with  the  wake  point  P^)  and  therefore  A ^  does  not 
require  any  additional  calculation. 

Now  the  wake  geometry  and  the  potential  discontinuity  are 
known  at  time  t*  +  dt  and  the  process  may  be  repeated. 


5 .2  Space  and  Time  Discretization:  First  Time  Step 

For  the  sake  of  clarity,  consider  first  the  case  of 
impulsive  start,  i.e.,  assume  that  for  t  <  0  the  rotor  was  at 
rest  in  a  fluid  also  at  rest.  Hence  at  time  t  *  0  there  is  no 
wake  surface  (the  potential  is  continuous  throughout  the  fluid 
region).  Equation  5.1  (without  wake  integral)  may  be  discretized 
by  dividing  the  surface  of  the  rotor  into  surface  panels  o- 
and  assuming  that  tjf  and  are  constant  within  each  panel.  J 

Nv 


2E(P,,0)  <f>(P*,0)  -  2  B.(P,.0) 

n*l  J  J 


Nb 

+  JLt  c.(p,,o)  u>. 

n=l  J  J 
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wh«re  and  ,  art  the  ▼aloat  of  ih  and  <j  on  the  j-th  panel  at 
tiae  t  *  0,  whereas  7 

bj(p...).JJ  i*, 

C,(P«.t)  -  ||  ?-<—-)  do  5.8 

J  JJ  _  a n  2nr 

By  imposing  the  condition  that  E q.  5.7  be  satisfied  at  the 
centroids  Pk  of  the  eleaenta  Oy.  one  obtains  (note  that  according 
to  Eq.  2.29  E(Pk)  »  1/2,  since  Pj.  is  a  regnlar  point  of  crk) 

"b  *b 

-  2  ^(0)  ‘f'j  +  £  <^(0)  fj  -  1 . J^)  5.9 

where 

B^U)  -  Bj  (Pk, t) 

<t )  -  Cj(Pk,t)  5.10 

It  shonld  be  eaphaaized  that  if  the  rotor  moves  with  rigid 
body  motion,  the  coefficients  B^,  and  Ckj  are  tiae  independent! 
Equation  5.9  is  a  system  of  ?fk'J  algebraic  equations  with  *Tt, 
unknowns  (the  values  of  :  are  known  from  the  boundary 
conditions).  J 

5.3 ,  lake  generation  Purina  First  Tiae  Step 

As  aentioned  above,  the  fraae  of  reference  is  assumed  to  be 
connected  with  the  air.  This  is  particularly  convenient  to 
discuss  the  wake  generation.  (The  actual  foraulation  used  in  the 
coaputer  program  is  for  a  frame  of  reference  connected  with  the 
blade:  the  two  are  related  through  a  rigid  body  rotation.  This 

point  is  discussed  later  in  this  section  where  it  is  indicated 
that  unstable  results  may  be  obtained  if  this  issue  is  not 
handled  carefully.) 

Between  time  t  •  0  and  time  t  «  At,  the  fluid  points  which 
were  at  the  trailing  edge  at  tiae  t  ■  0  move  into  a  new  position, 
which  (in  our  frame  of  reference,  i.e.,  a  frame  of  reference 
connected  with  the  undisturbed  air)  is  given  by  (see  Eq.  5.6) 

rAt 

Pw(At)  -  Pte(0)  +  9(Pw(t).t)  dt  5.11 


I 

An  analytic  expression  for  the  coefficients  is  siven  in  ~ef.  25. 
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From  a  discrete  (i.e.,  finite  element)  point  of  view,  the  new 
geometry  is  defined  by  the  location  of  the  wake  point  emanating 
from  the  trailing  edge  nodes  as  shown  in  Figure  5.1.  The 
location  of  these  points  may  be  is  obtained  by  approximating  the 
above  equation  as 

Pw(At)  *  Pte(0)  +  v(Pw(0).0)  At  5.12 

However,  the  calculation  of  v(Pv(0),0)  cannot  be  based  on  the  use 
of  Eq.  5.5.  The  reason  being  that  because  of  the  time 
discretization,  the  motion  history  is  actually  replaced  by  a 
sequence  of  sudden  starts  (see  Section  4.1).  This  implies  that 
in  the  time-discretized  formulation  a  discrete  vortex  is  formed 
at  each  time  step  and  therefore  the  velocity  as  given  by  Eq.  5.5 
goes  to  infinity  at  the  trailing  edge.  In  order  to  circumvent 
this  problem,  we  made  use  of  the  results  of  Section  4.4  and 
impose  that  the  relative  velocity  is  tangent  either  to  the  upper 
or  the  lower  side  of  the  rotor  blade  depending  upon  the  sign  of 
the  expression  in  Eqs.  4.12  and  4.14.  Then  we  used  standard 
finite-difference  expressions  to  determine  the  components  of  the 
perturbation  velocity  along  the  trailing  edge  (central- 
difference)  and  over  the  appropriate  side  of  the  blade  (e.g.,  the 
two-point  backward  difference  approximation  of  df^/dx  if  the  wake 
is  tangent  to  the  upper  side  of  the  blade). 

It  should  be  noted  that  as  mentioned  above,  the  new 
locations  of  the  wake  points  are  within  a  frame  of  reference 
connected  with  the  undisturbed  air.  Therefore  while  the  wake 
points  move,  the  blade  also  moves  into  its  new  position.  Hence  a 
row  of  wake  elements,  oQ  ,  has  been  generated  (see  Figure  5.1): 
these  elements  have  two  corners  on  the  trailing-edge  points  A  and 
B  whereas  the  other  two  corners  coincide  with  the  location  at 
time  t  *  At  of  the  fluid  points  that  coincided  with  the  blade 
points  A  and  B  at  time  t  =  0.  Bence,  at  time  t  *  At,  we  have  a 
row  of  wake  elements:  the  value  of  A(fa  assigned  to  the  elements 
a  is  the  difference  of  the  values  (evaluated  at  time  t  =  0)  of 
the  potential  if  at  the  centers  of  (the  upper  and  lower  blade 
elements  that  are  in  contact  with  o_  . 


5  .4 .  Generic  Time  Step 

Now,  at  time  t  3  At,  the  wake  geometry  and  the  values  of  A <f 
on  the  wake  are  known.  The  values  of  vf,  are  also  known.  Hence, 
Eq.  5.1,  discretized  as  (see  Eq.  5.7) 

Nb  Nb  N, 

£  <Hj  -  Ckj>  <fj  Ejjtj  +2  h.Afn  s-13 

may  be  used  to  evaluate  at  time  At.  In  Eq.  5.13,  B^j  and 
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Figure  5.1.  Wake  Generation  Between  Tiaes  t«  »  0  and  ti  ■  At 

i 


an  |Was  by  Eq.  5.10  whereas 


Fkn<‘> 


-ff  i-(i— )do| 

9a  Ip*-Pk 


5.14 


is  a  constant-doublet  integral  over  the  wake  element  on. 


Oaee  the  valaes  of  <f j  are  known.  the  saae  procedure 
described  ia  Sectioa  5.3  nay  be  repeated  to  calculate  the 
geoaetry  of  aew  wake  eleaeats  (aad  corresponding  values  for  Apn) 
geaeratsd  between  dt  aad  2dt.  Ia  addition,  the  new  location  of 
the  old  row  of  eleaeats  (those  generated  between  0  and  dt)  is 
obtained  as  follows:  evaluate  the  velocity  at  the  wake  points 
which  are  not  on  the  trailing  edge  using  Eq.  5.5  which  is 
discretized  as 


5.15 


where  q  spans  over  all  the  nodes  of  the  wake  surface  which  are 
not  on  the  trailing  edge,  whereas 


4« 

JJ 

•j 

~~i  do 
r  ! 

Ip*  *  \ 

?qj 

-1 

■  ■  — 

4  n 

K 

an  r 

dol 

Ip*  -  p, 

"SJ 

l, 

9  ,r 
— ( ~» ) 
3n  r 

doj 

!*•  -  P, 

5.16 


Then  calculate  the  new  locations  as 


Pq(t  +  dt)  -  Pq(t)  +  vq(t)dt  5.17 

Now  all  the  nodes  of  the  wake  surface  are  known.  Note  that  if 
the  node  nnaiberiag  of  the  wake  eleaeats  is  not  changed  frost  tiae 
dt  to  tiae  2dt  then,  according  to  Eq.  2.34 


d(f&  *  constant  in  tisM  5.18 

Hence.  the  new  wake  geoaetry  and  the  corresponding  values 
for  dfs  are  known  at  tiae  2dt  (as  aentioned  above,  the  geoaetry 
and  the  values  for  £tf  for  the  row  of  eleaeats  generated  between 
dt  aad  2dt  is  deterniaed  using  the  procedure  presented  in  Section 
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S.3)  and  the  process  may  be  repeated:  the  same  procedure  used  in 
the  time  step  between  At  and  2At  may  be  used  for  the  time  step 
between  tn  =>  nAt  and  tn+1  *  tn  +  At. 


3 .5 .  Maumricil  Stability 

Issues  related  to  numerical  instability  are  discussed  in 
this  section.  A  large  portion  of  this  research  effort  was 
completely  devoted  to  this  issue:  after  the  formulation  had  been 
completed  and  incorporated  into  a  computer  program,  the  results 
consistently  indicated  the  presence  of  numerical  instability.  It 
was  not  clear  at  that  time  whether  the  instability  correctly 
represented  a  physical  instability  (experimental  results  indicate 
that  the  wake  of  a  helicopter  rotor  in  hover  oscillates  and 
eventually  the  tip  vortex  breaks  down),  or  was  caused  by  the 
numerical  algorithm  (or  even  more  trivial,  a  bug  in  the  program). 
In  order  to  examine  the  possibility  of  a  bug  in  the  program,  a 
second  code  was  written  using  a  different  frame  of  reference, 
i.e.,  a  frame  of  reference  connected  with  the  undisturbed  air  (in 
the  first  code  the  frame  of  reference  is  connected  with  the 
rotor).  In  spite  of  the  equivalence  of  the  two  formulations,  the 
two  codes  gave  different  results.  The  source  for  the  difference 
was  traced  to  the  fact  that  in  the  first  code  the  new  location  of 
the  wake  point  was  evaluated  as 

Pw(t  +  At)  *  Pw(t)  +  (V*  +  vw)At  5.19 

where  vw  is  the  perturbation  velocity  whereas  V0  is  the  velocity 
of  the  undisturbed  fluid 


V0  *  2  i  t 


Unfortunately, 


APW  ■  Q  x  r  At  5 .20 

does  not  correspond  to  the  correct  rigid-body  displacement  (see 
Figure  5.3).  This  error  (which  goes  to  zero  as  At  goes  to  zero) 
seems  to  have  been  the  origin  of  the  instability:  the  results 
obtained  with  the  second  code  have  been  consistently  stable.  For 
the  hover  case,  it  is  actually  convenient  to  use  a  frame  of 
reference  connected  with  the  blade.  The  original  code  was 
therefore  modified  as  follows:  the  displacement  caused  by  the 
perturbation  velocity  is  calculated  then  the  whole^  wake  geometry 
is  rigidly  rotated  (around  a  vertical  azis,  i.e.,  fl  *  ftk)  as 

X  «  X,  cos  A0  -  y0  sin  AO 

y  *  Xo  sin  A9  +  y0  cos  A6 

z  -  z0  5.21 
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where  A6  *  CAt  whereas,  x0,  yo  ,  z«  are  the  coordinates  before 
rotation. 


5.6.  Finite-Core  Vortex 

Before  the  cause  of  the  instability  had  been  traced  to  the 
issue  discussed  in  Section  S.5,  the  concept  of  finite-core  vortex 
was  introduced  in  an  attempt  to  eliminate  the  instability.  (This 
concept  is  often  referred  to  as  artificial  viscosity:  we  believe 
that  this  name  is  misleading  and  prefer  the  name  finite-core 
vortex. ) 


The  essence  of  the  idea  is  that  concentrated  vortices  are 
introduced  in  the  process  of  spatial  discretization  of  the 
formulation.  The  integral  formulation  describes  the  wake  as  a 
doublet  layer  which  is  fully  equivalent  to  a  vortex  layer.  When 
the  wake  integral  is  divided  into  quadrilateral  panels  and  A<p  is 
assumed  to  be  constant  within  the  element,  we  are  in  effect 
replacing  a  continuous  distribution  of  vorticity  with 
concentrated  vortices  located  at  the  edges  of  the  panels.  Such 
an  approximation  therefore  introduces  infinite  induced  velocities 
in  the  vicinity  of  the  panels  edges.  If  the  wake  spirals  come 
too  close  to  each  other,  this  may  be  a  cause  for  instability. 

In  order  to  avoid  this  problem,  the  velocity  distribution 
induced  by  constant-doublet  panels  (i.e.,  concentrated  vortices) 
is  replaced  with  those  induced  by  a  finite-core  vortex.  Note 
that  we  considered  this  only  as  a  convenient  expedient  to 
eliminate  instabilities,  not  an  effort  to  improve  the  model  by 
making  it  closer  to  physical  reality  (in  order  to  get  closer  to 
physical  reality  we  would  rather  introduce  a  higher-order  finite- 
element  representation  for  the  distribution  of  the  potential 
discontinuity  in  place  of  the  constant-doublet  approximation). 

In  order  to  accomplish  that,  the  following  scheme  was  used 
in  our  computer  code.  The  velocity  induced  by  the  edge  PXP,  of  a 
constant  doublet  panel  is 


v 


1  qixq» 

2n  lqixq*|* 


qx • qi~qi  •  q* 
( - 


Iqxl 


q*  * q*-qi • q» 
- ) 

Iqi  I 


where  q^ 


-  P*  (see  Figure  5.4). 


5.22 
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This  is  exactly  equal  to  the  expression  obtained  froa  the  Biot 
and  Savart  law  for  the  Telocity  'induced*  by  a  vortex  line. 
Introducing  a  fraae  of  reference  such  that  the  vortex-line 
coincides  with  the  z-axis  and  that  the  x-axis  goes  through  P*, 
Eq.  5.22  yields  (see  Figure  5.2) 

▼  - —  lx»(x  +ii  )  '  +z»(x  +  z»)  '  ]j  5.23 


As  P*  approaches  the  vortex  line  (i.e.,  as  x  tends  to  zero),  this 
expression  yields 


v 


1  - 
-  j 

2nx 
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(which  is  the  velocity  induced  by  an  infinite  vortex)  if  P* 
approaches  a  point  between  Pa  and  Pa  (i.e..  zx  <  0  and  i,  >  0) 
and 


v  -  0  5.25 

If  P*  approaches  a  point  of  the  line  PiPs  outside  the  segaent 
paPa  ( i.e.,  if  Zi  <  fa  <  0  or  za  >  za  >  0).  An  infinite  vortex 
with  a  finite  core  of  radius  e  yields  a  velocity 


1  - 
a  -  j 

2irx 
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.  where 

a  -  1  x  >  s 

“  x*/e*  x  <  s  5.27 

Hence,  going  bach  to  vector  notations.  Equation  5.26  say  be 
rewritten  as 


1  5i*qa 
R  lq.xi.1* 


qi*qx-qa*qa 


qa*qs-qx*qa 
- ) 

Us  I 
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if  x  >  s  and 


1  1  5a.*5s 


qx’qz-qx’qs 
( - 


Iqal 


qa’qa-qa’qa 

- ) 

Iq.l 
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if  x  <  s,  where 


I  'Ll  x  q*l 

x  -  — -  5.30 

I  ii  -  5a  I 
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For  the  take  of  clarity,  the  formulation  has  been  presented 
for  the  ease  of  start  from  rest,  that  is,  for  tine  t  <  0,  the 
rotor  has  been  as  sued  to  be  at  rest  and  surrounded  by  fluid  also 
at  rest. 

It  should  be  apparent  at  this  point  that  such  an  assuption 
is  not  essential  to  the  formulation.  The  method  of  solution 
requires  only  that  the  geometry  wake  surface  and  the  potential- 
discontinuity  distribution  oyer  the  surface  of  the  wake  to  be 
known  at  time  t  Then,  the  potential  distribution  may  be 
evaluated  using  Eq.  S.13.  Once  the  potential  distribution  is 
known,  Eq.  5.15  may  be  used  to  calculate  the  velocity  of  the 
nodes  of  the  wake  and  hence  their  new  location  (the  row  of  the 
new  wake  elements  generated  between  tQ  and  tQ  +  At  is  obtained 
with  the  procedure  given  in  Section  5.3). 

As  mentioned  above,  the  frame  of  reference  connected  with 
the  undisturbed  air  is  more  convenient  from  a  conceptual  point  of 
view  even  though  the  calculations  are  more  conveniently  presented 
in  a  frame  of  reference  connected  with  the  rotor.  Also,  the 
issue  of  artificial  viscosity  discussed  in  Section  5.6  does  not 
appear  to  be  essential  to  the  formulation:  small  values  of  e  give 
the  same  results  as  s*0. 

The  last  major  issue  to  be  discussed  is  that  of  the  wake 
truncation:  as  the  number  of  time  steps  grows,  the  length  of  the 
wake  also  grows.  This  implies  that  the  computer  time  per  time 
step  also  grows.  In  order  to  keep  computer  time  within 
reasonable  bounds,  it  is  necessary  to  obtain  a  simplified  model 
for  the  remote  element  of  the  wake.  Vhile  sophisticated 
intermediate-and  far-wake  models  have  been  introduced  for  the 
hover  case  (Refs.  40  and  41),  these  models  require  ad-hoc 
assumptions  based  on  empirical  data.  Since  the  objective  of  the 
present  work  is  to  develop  a  method  which  may  be  used  to  study 
problems  for  which  such  data  does  not  exist  (such  as 
maneuvering),  it  would  have  been  inappropriate  to  introduce  any 
of  the  above  far-wake  models  or,  for  that  matter,  any  model  based 
on  experimental  data.  For  this  reason  in  the  results  presented 
here,  the  wake  is  simply  truncated  after  a  certain  number  of 
spirals.  The  implication  of  this  procedure  is  that  the  last 
few  spirals  are  to  be  considered  as  modelling  of  the  far  wake 
effects.  As  indicated  by  the  results  presented  in  Section  6,  this 
is  an  expenaive  approach  to  the  problem  (the  case  presented  in 
Section  6  requires  approximately  eight  hours  of  CPU  time  on  an 
IBM  370/168).  There  is  need  to  develop  a  less  expensive  approsch 
to  the  far  wake  modelling.  However,  such  a  model  should  be  based 
on  first  principle  rather  than  empirical  data,  if  the  methodology 
proposed  here  is  to  be  used  independently  of  the  experimental 
analysis. 
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SECTION  6 


NUMERICAL  RESULTS 


In  order  to  validate  the  theory  presented  above,  the 
numerical  algorithm  was  implemented  in  a  computer  program.  The 
results  obtained  with  this  program  and  the  comparison  with 
existing  data  are  presented  in  this  section. 


6.1,  Choice  of  the  Test  Case 

It  is  a  generally  accepted  opinion  in  helicopter-rotor 
aerodynamics  that  the  wake  roll-up  problem  is  harder  to  solve  for 
hover  than  for  forward  flight:  this  is  because  the  wake  spirals 
are  closer  to  each  other  in  the  hover  case.  Also  the  hover  case 
seems  to  be  the  only  one  for  which  satisfactory  results  exist. 
Therefore,  in  order  to  test  our  formulation  we  have  started  by 
studying  a  hover  case.  In  order  to  validate  the  time-domain 
algorithm  the  hover  case  was  studied  through  a  time-accurate 
transient  response  analysis.  These  steady  state  results  are  the 
only  ones  presented  here.  We  believe  that  the  validation  of  the 
formulation  will  be  satisfactory  only  if  more  extensive  results, 
(including  forward  flight  results  now  under  consideration)  will 
confirm  the  results  presented  here. 

In  particular,  we  chose  the  case  studied  by  Rao  and  Schatzle 
(Ref.  7)  for  several  reasons,  the  most  important  of  which  is  that 
their  formulation  (lifting  surface  with  prescribed  wake)  is  based 
on  first  principles  (no  ad-hoc  assumption  is  used  except  for  the 
wake  geometry  and  zero-thickness  blade)  and  yields  results  which 
are  in  excellent  agreement  with  the  experimental  ones  of  Bartsch 
(Ref.  37).  Also  important  in  Ref. 7  is  the  comparison  between 
classical-  and  Landgrebe-wake  analysis:  the  latter  shows  a 
marked  spike  near  the  tip  in  the  section  lift  distribution,  in 
full  agreement  with  experimental  results  (the  spike  is  absent  in 
the  classical  wake  results).  These  results  not  only  demonstrate 
the  importance  of  the  wake  geometry  for  the  prediction  of  the 
section  lift  distribution,  but  also  are  important  in  the 
interpretation  of  our  convergence  scheme. 


6.2.  Description  of  Test  Case 

As  mentioned  above,  all  the  results  presented  in  this 
report  for  a  case  considered  by  Rao  and  Schatzle:  an  isolated 
rotor,  with  tip  radius  Rj  »  17.5',  cut-out  radius  YCQ  =  2.33', 
chord  c  =  1 .083',  collective  pitch  angle  0  =  10.61  and  twist 

angle  0X  *  -5  .  The  a-  ular  speed  is  £1  =  355  r.p.n. 

For  all  the  results,  the  initial  wake  geometry  is  a 
classical  wake 
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x  =  *0  cos  0  -  y0  sia  0 


y  3  y0  sin  0  +  x0  sin  0 
z  =  z0  +  kR0 

where  k  =  J  Cj./2  and  Cj.  =  0,00186  (this  is  the  value  obtained  by 
Rao  and  Schatzle). 

All  the  results  were  obtained  using  three  elements  in  the 
chord  directions  and  seven  in  the  span  directions  for  a  total  of 
twenty-one  elements  on  each  side  of  the  blade.  A  convergence 
analysis  presented  in  Ref.  32  indicates  that  this  is  sufficient 
to  obtain  relatively  converged  results.  For  the  classical  wake 
analysis  the  time  step  is  At  3  T/12  where  T  is  the  period:  this 
yields  twelve  elements  in  the  ‘circumferential*  direction  per 
each  wake  spiral  (there  are  seven  elements  in  the  radial 
directions  because  that  is  the  number  of  elements  on  the  blade  in 
the  spanwise  direction). 

These  data  are  summarized  in  Table  A. 


6.3  Numerical  Results 

In  order  to  illustrate  the  type  of  problems  encountered  in 
the  effort  presented  here,  consider  Figures  6.1. a  to  6.1.p  which 
depicts  the  vertical  displacement,  z,  as  a  function  of  the 
azimuth  angle  0X  (for  different  time  steps,  more  precisely  for 
tQ/At  3  15  and  21  through  35)  for  the  last  vortex  line,  i.e.,  the 
vortex-line  emanating  from  the  tip  point  of  the  trailing  edge  (as 
clarified  later  it  would  be  misleading  to  refer  to  this  vortex  as 
the  tip  vortex). 

It  may  be  noted  that  at  time-step  15  the  vortex  line  appears 
to  have  the  expected  behavior.  However  at  time  step  20,  the 
vortex  line  shows  a  hump  which  becomes  more  marked  as  in  the 
following  time  steps  but  it  seems  to  be  pushed  downward  in  time. 
The  phenomenon  is  more  clearly  illustrated  in  Figure  6. 2. a  to 
6.2.d  which  depicts  the  vertical  displacement,  z,  as  a  function 
of  the  radial  position  (for  different  time  steps,  tn/At  3  15,  22, 
23,  25)  also  for  the  last  vortex  line.  It  is  apparent  that  the 
last  spirals  tend  to  move  outward.  This  is  caused  by  the  fact 
that  the  wake  is  truncated  (the  ‘following’  spiral  would  have  a 
'restraining'  effect  on  the  last  spiral:  in  its  absence  the  last 
spiral  tends  to  move  outward,  this  point  is  clearer  after 
discussion  of  Figure  6,5). 
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Single  Sot  or  Blade 

Geometry  and  Flight  Conditions  of  Rotor: 

Number  of  Spirals  5  or  7  at  indicated  in  captions 
Number  of  Elements  on  Blade  (3  x  7) (one  side) 
Number  of  Elements  per  Spiral  (12) 

R  -  17.5  ft.  (span) 

o  *  1.083  ft.  (chord) 

Tco  ■  2.33  ft.  (root  cut-out) 

®r  *  10.61*  (collective  pitch  angle) 

*  -5*  (twist  angle) 

Ct  *  0.00186  (thrust  coefficient) 
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ling  edge  at  tlaeetep  H  33(for  caaa  defined  in  Table  A). 
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As  Mentioned  in  Section  5.8,  the  reason  no  far-wake  aodel 
was  introduced  is  beoanse  sack  a  aodel  does  not  exist  for 
arbitrary  notion  (although  the  far-wake  aodel  introdnced  in  Ref. 
40  and  41  for  the  hover  case  conld  have  been  nsed  here  in  order 
to  iaprove  these  specific  results). 

In  spite  of  the  fact  that  the  wake  geoaetry  does  not  seen  to 
have  reached  convergence,  the  lift  distribution  appears  to  be 
converged:  no  appreciable  changes  occur  between  tiae  steps  40 
and  SO.  This  is  shown  in  Figures  6.3. a  and  6.3. b  for  the  five- 
and  seven-spiral  models  respectively:  the  converged  five-  and 
seven-spiral  results  are  coapared  in  Figure  6.4.  Ve  believe  that 
the  explanation  for  this  seeaingly  inconsistent  behavior  can  be 
understood  froa  the  results  shown  in  Figure  6.5  which  shows  the 
vertical  displaceaent  of  the  last  vortex  line  as  a  function  of 
the  axiauth  angle  at  the  last  tiae  step  considered  here  (that  is, 
at  tiae  step  50):  the  three  lines  correspond  to  the  three-, 
five-,  and  seven-spiral  aodels  reapectively.  It  aay  be  seen  that 
the  first  two  spirals  are  very  close  for  all  three  cases  (whereas 
the  first  three  spirals  are  in  good  agreeaent  for  the  five-  and 
seven-spiral  aodels).  Since  the  first  one  or  two  spirals  have 
the  strongest  inpact  on  the  section  lift  distribution,  it  is  no 
surprise  that  the  five-  and  seven-spiral  aodels  are  in  good 
agreeaent  on  the  section-lift  distribution. 

Finally,  the  results  obtained  with  the  seven-spiral  wake  are 
coapared  against  existing  data  in  Figures  6.6,  6.7,  and  6.8. 
Figure  6.6  shows  a  cross  section  of  the  wake  (first  two  spirals 
only)  at  90*  behind  the  trailing  edge.  Also  shown  in  Figure  6.6 
are  the  location  of  the  tip  vortex  and  of  the  vortex  sheet  as 
predicted  by  Landgrebe's  generalized  wake  aodel.  Note  that 
Landgrebe's  aodel  coaes  froa  the  experiaental  data  and  therefore 
the  tip  vortex  is  not  necessarily  the  location  of  the  last 
vortex,  but  rather  the  'center  of  aass’  of  the  rolled-up  portion 
of  the  vortex  sheet.  Taking  this  into  account,  we  do  consider 
this  coaparison  to  be  very  satisfactory  especially  if  one 
considers  the  low  nuaber  of  eleaents  used  to  describe  the  blade 
and  its  wake:  auoh  stronger  roll-up  is  expected  if  a  higher 
nuaber  of  eleaents  is  used.  (It  aay  be  worth  noting  that 
Landgrebe  aodel  is  only  an  approxiaate  interpolation  of  the 
experiaental  data.)  Siailar  good  agreeaents  are  shown  in  Figure 
6.7  in  which  the  radial  location  of  the  last  vortex  as  a  function 
of  the  aziauth  0  is  coapared  to  the  radial  location  of  the  tip 
vortex  in  Landgrebe's  aodel. 

Finally,  Figure  6.8  shows  a  coaparison  of  our  results  with 
those  of  Rao  and  Schatzle  (Ref.  7).  As  Mentioned  in  Section  6.1, 
their  results  for  a  four-bladed  rotor  are  in  excellent  agreeaent 
with  the  experiaental  results  of  Bartsch  (Ref.  37).  Again,  we 
consider  that  the  agreeaent  is  satisfactory  if  one  considers  the 
low  nuaber  of  eleaents  used  in  the  analysis  and  that  Rao  and 
Schatzle  results  are  obtained  with  a  prescribed  wake.  (It  aay  be 
worth  noting  that  our  results  are  in  excellent  agreeaent  with 
their  results  for  classical-wake  analysis,  see  Ref.  32.) 


Figure  6. 3. a.  Convergence  analysis:  spanwise  sectional  lift  distribution, 
(lb/ft),  as  a  function  of  radial  distance,  r/R,  for  five  spiral  wake  (for 
cite  defined  in  Table  A) • 
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M,  Cgmsnjjj  Numerical  Stability 

From  the  numerical  reanlta  presented  in  this  Section  one  may 
conclude  that  the  algorithm  ia  capable  of  reproducing  the  correct 
trend  in  wake  rollup  and  pressure  distribution.  The  discrepancy 
between  our  results  and  the  eziating  ones  may  be  due  to  either 
the  physical  approximation  (i.e.,  inviseid  flow)  or  numerical 
approximation.  An  analysis  of  convergence  is  needed  in  order  to 
discriminate  between  the  two.  (The  irregular  behaviour  it  0  * 
60  is  probably  due  to  the  trailing  edge  condition  discussed  in 
Section  5J  and  deaerves  further  attention). 

The  main  accomplishment  however,  is  that  the  numerical 
results  indicate  that  the  algorithm  appears  to  be  free  of 
numerical  instabilities,  even  though  no  ad-hoc  assumption  (such 
as  prescribed  radial  contraction)  has  been  used. 

More  precisely  we  believe  that  the  instability  in  the  last 
few  spirals  on  the  wake  is  due  to  the  truncation  of  the  wake  and 
should  not  be  thought  of  as  a  numerical  instability  in  the 
classical  sense:  in  such  a  case  the  vortex  line  would  depart  from 
a  amooth-behavior  spiral  with  a  disturbance  that  oscillates  and 
grows  in  space  and  time  such  as  that  reported  by  Summa  (Ref.  6). 
All  of  our  results  are  very  smooth:  an  illuatrative  example  of 
such  smooth  behavior  is  presented  in  Figure  6.9,  which  shows  the 
geometry  of  the  vortex-line  emanating  from  the  tip  point  of  the 
trailing  edge  for  all  the  spirals  (for  the  five-spiral  case)  at 
time  step  n  =  45. 

Ve  believe  that  this  is  the  first  time  that  such  an 
accomplishment  has  been  reported. 
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point  of  the  trailing  edge  for  all  the  spirals  (for  the  five- 
spiral  case)  at  time  step,  n- 45  (for  case  defined  in  Table  A). 


SECTION  7 


concludin':  renarss 


A  general  formulation  for  the  aerodynamic  analysis  of  the 
unsteady  incompressible  potential  flow  around  a  helicopter  rotor 
has  been  present ed  (the  extension  to  compressible  flows  is 
outlined  in  Appendix  A).  New  theoretical  results  incorporated  in 
the  formulation  include  a  thorough  discussion  of  the  wake 
dynamics  (Section  2)  and  of  the  wake  generation  (Section  4),  both 
examined  from  very  basic  principles.  In  particular  the  trailing 
edge  condition  is  examined  in  detail.  It  is  shown  that  ITutta 
condition  (that  the  vorticity  distribution  does  not  go  to 
infinity  at  the  trailing  edge),  necessary  to  insure  uniqueness  of 
the  solution  for  two  dimensional  flows,  is  not  needed  for 
three  dimensional  flows  (in  simply  connected  regions).  On  the 
other  band  Joukowski  hypothesis  (that  concentrated  vortices,  if 
formed,  are  immediately  shed  off  the  trailing  edge)  is  used  in 
the  discussion  of  the  wake  generation.  The  numerical  algorithm 
is  also  presented:  it  is  shows  how  the  tine-discretization  of  the 
problem  yields  a  series  of  'sudden  starts'  with  the  consequence 
that  discrete  vortices  are  generated  (and  shed)  at  the  trailing 
edge. 


Numerical  results  are  then  introduced.  The  main  conclusions 
which  may  be  drawn  from  these  results  are  that 

1.  the  algorithm  appears  to  be  numerically  stable 

2.  seven  spiral  wakes  are  needed  to  obtain  ’good  behavior’ 
for  the  first  two  or  three  spirals 

3.  the  finite-core  assumption  (artificial  viscosity)  does 
not  seem  to  affect  the  results 

4.  the  wake-geometry  results  are  in  good  agreement  with  the 
empirical  analytical  model  of  Landgrebe  (Ref.  33) 

5.  the  section  lift  distribution  is  in  good  agreement  with 
the  results  of  Rao  and  Schatzle  (Ref.  7) 

The  most  important  new  development  is  that  to  our  knowledge 
this  is  the  first  time  that  results  have  been  obtained  without 
the  need  for  an  empirical  input  to  avoid  numerical  instability. 

Although  the  validation  has  been  obtained  only  for  a  rotor 
in  hover,  the  formulation  is  quite  general  (the  main  limitations 
being  irrotationab ility  and  incompressibility)  and  applicable,  in 
particular,  to  a  rotor  in  forward  flight. 
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Additional  work  it  recoaaended  ia  the  following  treat: 


1.  Convergence  analvtit :  it  it  expected  that  a  ttronger 
roll-up  would  be  obtained  by  uting  a  larger  number  of 
eleaenta  in  tbe  radial  direction  (tbit  in  turn  would 
affect  the  aection-lift  diatribution). 

2.  !»is  truncation:  it  it  recoaaended  that  toae 

intermediate-  and  far-wake  model  be  introduced  for 
the  purpote  of  reducing  the  nuaber  of  wake  tpiralt 
(and  hence  the  CPC  tiae).  However,  at  mentioned  in 
Section  S.7,  tuch  aodelt  thould  be  bated  on  firat 
principlet  rather  than  on  eapirical  data,  if  the 
objective  of  the  methodology  it  to  uae  it  for  the 
calculation  of  generalized  waket. 

3.  Validation:  continue  the  validation  of  the  formulation 
by  applying  it  to  additional  hover  catet  and  then  to 
forward  flight  catet  (it  ia  hoped  that  comprehent ive 
experimental  reaulta  for  forward-flight  become 
available  in  the  near  future). 
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APPENDIX  A 


FORMULATION  FOR  COMPRESSIBLE  FLOWS 

Tie  integral  formulation  of  Section  3  is  extended  here  to 
the  case  of  compressible  flows.  The  frame  of  reference  is 
assumed  to  have  arbitrary  motion.  The  surface  is  assumed  to  be 
moving  with  respect  to  the  frame  of  reference  in  order  to 
accommodate  structural  deformations  as  well  as  wake  roll-up. 
However,  for  the  sake  of  simplicity,  such  motion  is  assumed  to  be 
small.  The  general  case  is  considered  in  Refs.  63  and  64.  The 
formulation  is  an  extension  of  that  introduced  in  Ref.  65  for 
acoustics . 


Equivalent  Problems 

The  equation  for  the  velocity  potential  in  a  frame  of 
reference  connected  with  the  undisturbed  air  is  given  by 
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A.l 


where  contains  all  nonlinear  terms.  The  boundary  conditions 
represent  the  zero  flow  at  infinity 

f  “  0  at  ***'  A. 2 

and  the  no- flow- through  condition  on  the  surface  boundary 


•  . 
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In  addition. 

D_ 

—  Let  -  0 

Dr  ^ 

on  wake 

In  order  to  simplify  the  derivation  of  the  Green's  theorem  for  a 
frame  of  reference  having  arbitrary  motion,  it  is  convenient  to 
extend  the  problem  to  the  whole  space  by  introducing  the  function 
^  *  E ij>  where,  for  a  surface  o  surrounding  both  body  and  wake 

E  *  1  inside  a 

■  0  outside  o  A. 5 

so  that 

f  -  f  inside  a 

-  0  outside  9  A.6 
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Tote  that  (ji  is  defined  in  the  whole  space  and  satisfies  the 
enua t ion 
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The  presence  of  the  7E  and  dE/dr  terns  introduce  source  layers 
which  act  only  on  the  surface  a  (which  is  not  to  be  considered  as 
a  boundary  of  the  domain  of  validity  of  the  equation  which  is  the 
infinite  space).  The  only  boundary  is  at  infinity  where  v;e 
specify  cp  =  0.  By  denoting  the  right  hand  of  Equation  A.7  by 
Eq.  A.7  nay  be  rewritten  as 


9  a 


3t 
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whereas,  according  to  Eq.  A. 2,  the  boundary  condition  at 
infinity,  is 

Cf>  =  0  at  <o  A. 9 


Equation  A. 3  subject  to  boundary  condition  A. 9  is  equivalent  to 
Equation  A.l  to  A. 4  in  the  sense  that  if  a  function  satisfies 
E  q  s .  A. 3  and  A. 9,  it  also  satisfies  E  q  s .  A.l  to  A. 4.  The 
solution  to  Equations  A.8  and  A. 9  is 


»  — 


A.  10 


where 


O(l.r) 


6 (x*-t-9) 


4np 


A.  11 


(with  p  =  and  9  =  p/a*)  is  the  well  known  Oreer.'s  Function 

for  the  wave  operator. 


Transfornat  ion  to  the  "oving  Trane  of  Ee  ferenc e 


"ere  we  introduce  a  coordinate  system  (x,t)  moving  in  rigid- 
body  notion  (e.g.,  connected  to  the  rotor  if  the  rotor  moves 
with  rigid  body  motion).  v.'e  nake  use  of  a  transformation  from 
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the  fixed  frame  of  reference  (?,t)  to  the  moving  frame  of 
reference  (x,t) 


dt 


V(x, t) 


A. 12 


$(x,0)  “  x 


Then 

5  ( x.  t )  *  x  + 
T  -  t 

From  Equation  (6)  we  have 
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X(x.t)  S(t-t,+6(i,t)) 
4rrp(x,t) 


dx  dt 
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A. 14 
A. 15 


A. 16 


where 

p(i.t)  -  l?(x,t)  -  5(x,>t*)l  A. 17 


(Since  the  velocity  V  corresponds  to  rigid  body  motion,  the 
transformation  is  a  length  preserving  and  hence  the  Jacobian  is 
equal  to  one.)  Recalling  a  result  from  theory  of  distributions 
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f(t)  8  [g(t)l  dt 
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(where  t^  are  the  roots  of  g(t)  ■  0)  and  integrating  with  respect 
to  the  time  variable,  one  obtains 


4 if  <f  (**,  t*) 

where  rfT(x,  t . )  „  p(x,ti)  !  g'ft  £) ! 
the  solutions  of 

g(t)  ■ 


(I, ti> / rn(x, ti> ]  dx  A. 19 

,  whereas  the  values  of  tj  are 
>  0  A. 20 


where 

g(t)  ■  t*  -  t  -  !  ?(x,t)  -  f(x*,t*)l /a®  A.21 

Ve  note  that  the  Eq.  A.19  is  the  desired  integral  representation 
for  the  velocity  potential  in  the  frame  of  reference  moving  with 
the  body. 
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Integral  donation  for  Non- 


r.laid  n.otor 


Mote  that  if  %*  C  (i.e.,  if  the  nonlinear  terns  in  the 
differential  equation  for  the  velocity  potential,  Eq.  A.l,  are 
negligible)  then  the  right  hand  side  of  Eq.  k.l ,  is  equal  to 
zero  except  on  the  surface  o  (since  E  is  constant  both  outside 
and  inside  o).  Ilence,  if  %m  0,  the  right  hand  side  of  Eq.  A.19 
can  be  rewritten  in  the  forn  of  the  surface  integrals. 

For  the  sake  of  clarity,  consider  first  the  case  of  a  non¬ 
lifting  rotor  in  rigid  body  notion.  Let  the  coordinate  system  ? 
coincide  with  the  Cartesian  coordinates  of  a  Cartesian  frame  of 
reference  rigidly  connected  with  the  rotor.  In  this  frame  of 
reference 


3E 

dt 
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In  addition,  assune  that  the  rotor  speed  is  always  subsonic  so 
that  Eq.  A. 20  has  only  one  root,  t»T. 


Note  that,  for  any  function  f 
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do 


and 


JJJ  d  f-ifVK)  dV  -  P*  'JJJ  <3f$E  dV 

-  Vjj  GfS  do 

Performing  the  calculations,  one  obtains  (see  P.ef. 
details) 
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In  Eq.  A.25,  (  ]T  indicates  evaluation  at  retarded  tine  t-T  (such 
that  g(t)  «  t  -  t,  +  !$(x,t)  -  ?(x*,t*) l/a.  -  0).  Equation  A.25 
is  the  desired  integral  representation.  In  the  limit,  as  T* 
approaches  the  surface  o,  one  obtains  an  integral  equation  for 
(see  Section  3).  The  numerical  solution  of  such  an  equation  Is 
similar  to  that  given  in  Section  5. 


lotof/F«nlt«>  CotflMMtiw 


Hast  eoaaidar  tk«  tin  of  rotor/ fact  lag*  eoaf  igaratioa  ia 
vhieh  botk  tki  rotor  tad  tki  faaalaga  mot*  ia  arbitrary  bat 
rigid-*ody  aotioaa.  Alao  for  siaplieity,  aaaaaa  that  the  viki 
reaaiaa  where  it  ia  geaerated:  thia  ia  a  raaaoaabla  aaaaaptioa 

whoa  (ia  a  fraao  of  roforaaea  eoaaoctad  with  tha  aadiatarbad  air) 
tha  walooity  of  tha  flaid  ia  aaall  eoaparad  to  that  of  tha  rotor 
faaalaga  aoaf igaratioaa  (thia  aaaaaptioa  ia  raaowad  ia  tha 
aaalyaia  of  Kaf.  «3).  Haaea  tha  aarfaea  e  eaa  ba  brokaa  into 
thraa  aarfaoaa:  tha  aarfaaa  of  tha  rotor.  ot,  tha  aarfaea  of  tha 
faaalaga  Of,  aad  tha  aarfaea  of  tha  waka,  ov.  For  aaeh  of  thaea 
aarfaeaa  thara  aziata  a  fraao  of  rafaraaea  vhieh  ia  rigidly 
eoaaeetad  with  tha  aarfaea. 

Nazt  eoaaidar  tha  volaae  oa  tha  right  head  aide  of  Eq.  A.18 
aad  eoaaidar  thraa  'thia'  volaaes  eoataiaiag  or,  af  aag 
raapaetiwaly  (tha  portioa  of  tha  iatagraad  eoataiaiag  derivative* 
of  B  ia  eqaal  to  aaro  ia  tha  raat  of  tha  volaae).  Nazt  assaae 
that  1  ia  aaeh  volaae  eoiaeidaa  with  the  fraao  of  refereaea 
eoaaaetad  with  tha  eorreepoadiag  aarfaea.  Followiag  the  aaaa 
proeadara  oatliaad  above  (bat  loeatiag  eaeh  aarfaea  iategral 
iadepeadeatly).  oaa  obtaiaa  aa  azpaaaioa  aiailar  to  Bq.  A.25  with 
aaeh  iategral  raplaead  with  tha  aaa  of  thraa  iatagrals  over 
•t,  Of  aad  ov  reapeetivaly. 

If  the  aotioa  of  tha  aarfaeaa  with  raapaet  to  'their  fraae 
of  rafaraaea'  i*  aaall,  Bq.  A.23  is  atill  valid  bat  aaeh  aotioa 
'above  ap'  ia  the  boaadary  eoaditioaa  for  fly/9a.  Tha  ease  of 
eoaplataly  arbitrary  aotioa  is  diseasaod  ia  Kaf.  63  (iaeladiag 
aoaliaaar  toras  aad  raaoviag  shocks):  tha  darivatioa  of  the 
eqaatioaa  ia  vary  ooaplieatad  bat  tha  fiaal  resalta  are  slightly 
aoro  eoaploz  thaa  the  ones  praaeatod  hare. 


